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ABSTRACT 



A pulse-echo system using an electrostatic transducer and 
steel reflector was operated over the frequency range 24-192 Kcps 
in the sea and in the laboratory. Excess attenuations from 0.15 to 
1.1 db/M due to absorption and scattering by inhomogeneities were 
made over the full frequency range with a maximum error of less than 
0.05 db/meter. Backscattering at sea was detected in the frequency 
range 24-60 Kcps with an average scattering cross-section per unit 
volume of 400 cm^/m^ at a depth of 10 feet. Scattering was not 
detected gbbve 60 Kcps. Variation of sound velocity with frequency 
was measured from 24-192 Kcps. The accuracy of these values was 0.03 
percent for frequencies over 40 Kcps. Results are compared and it is 
postulated that; 

(1) below 60 Kcps the absorption and scattering was due to 
microbubbles of radius greater than 50 microns, and 

(2) above 60 Kcps minute plankton caused absorption without 
significant scattering. 
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1. Introduction 



It has been found that the behavior of sound upon being 
reflected from the surface of the ocean does not agree with theory, 
particularly for angles of incidence away from the normal. One 
of the possible explanations for this disagreement is the presence 
of small bubbles or other particulate matter near the surface. 

Theory and laboratory experience have shown that a small concentra- 
tion of bubbles in water would be particularly effective in causing 
backscattering, attenuation, and changes in sound velocity. It is 
the purpose of this project to experimentally determine the presence 
of microbubbles and other scatterers and absorbers near the ocean 
surface. 

A basic pulse-echo system [28^] consisting of an open arrangement 
of an electrostatic piston transducer facing a reflector was used 
to determine attenuation, scattering and velocity dispersion at two 
field locations described in Appendix VI. The three types of 
measurements were then used in order to interpret the sources of 
absorption and scattering of sound at frequencies ranging from 
24 Kcps to 200 Kcps, in sea states one and two, and primarily at 
a depth of ten feet. 
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2 . Instrumentati on 



The instrumentation used in this project consists of four 
main sections; a transmitting section, a receiving section, a 
transmit-receive switching section, and a transducer section. 

Each will be described briefly. 

The transmit section takes a sinusoidal signal, gates it to 
form wave pulses, and then amplifies these wave pulses. The wave 
pulses are added to a direct current polarizing voltage and then 
applied to the electrostatic transducer through the transmit- 
receive switching section. Typical transmitter outputs are shown 
in figure 2-1. A block diagram of the transmit section may be 
found in figure 2-2. 

The transmit-receive switching section utilizes a signal diode 
network and an electronic switch to isolate the transmitter from 
the reciever. Figure 2-3 is a block diagram of this section. 

Figure 2-4 shows the receiver output while transmitting . 

The receiving section consists of an electronic filter and 
oscilloscope display with simplification before and after the filter. 
Figure 2-5 is an example of a receiver output display. Figure 2^6 
is a block diagram of the receiving system. 

The transducer section consists of a 24 cm diameter wide-band 
electrostatic piston transducer and a parallel 30 cm steel 
reflector located 76 centimeters along the axis from the transducer. 

A complete description of the design, construction and characteristics 
of the electrostatic transducer is given in Appendix I. 

Figure 2-7 shows the electronics set up at the Naval Electronics 
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Laboratory Tower. Manufacturer's specifications for the commercial 
equipment used in the entire system may be found in Appendix V. 

The overall system characteristics are given in Appendix IV. 
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Figure 2-2 
Transmit Section. 
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Figure 2-3 

Transmit-Receive Switching Section 
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Figure 2-6 
Receive Section 
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3 . Attenuation Measurements 



a. Basic Principle , 

An Electrostatic transducer was gated to emit a short 
pulse and then electronically switched to receive the successive 
echoes returned from a reflector placed 76.4 cms. in front of it. 
Photographs of the pulse echo return are shown in Fig. 3.1. In all 
cases the gain of the bscilloscope was adjusted so that the 
amplitude of the second echo was approximately 10 cm. when measured 
on the oscilloscope face. 

Successive echoes were reduced in amplitude by three types 
of loss; divergence, imperfect reflection at the reflector and the 
transducer and attenuation in the medium. This attenuation in the 
medium was due to scattering and absorption by inhomogeneities, in 
addition to the conventional absorption of the medium itself. The 
object was to measure the excess attenuation loss in the sea. For 
all media studied, in the frequency range 24-192 Kcps, this medium 
loss was smaller than the reference system loss due to reflection 
and divergence. 

At the lower frequencies, below approximately 60 Kcps, 
divergence loss was very great so that the number of pulses which 
could be measured was smaller than that desirable for graphical work. 
However graphs were drawn for frequencies as low as 24 Kcps and used 
for the calculation of attenuation. 

Fig. 3.1 shows photographs from the minimum usable 
frequency of 24 Kcps to the maximum frequency at which it was 
possible to operate the equipment satisfactorily, 240 Kcps. 
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c) 96K cps. 



d) ^40 Kcps. 



Fig. 



Photographs obtained during the calibra tion of the system . 
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Higher frequency operation will require the use of an amplifier 
capable of delivering an undistorted input to the transducer at 
the frequency of interest and a filter capable of operation at that 
frequency. As can be seen from Fig. 3.1 d, noise made the accurate 
measurement of pulse amplitude increasingly difficult above 200 
Kcps. 

b. Interference Effects 

Interference effects are illustrated by Fig. 3.2 which 
shows pulse amplitude (logarithmic scale) against echo number 
(linear scale). These effects are superimposed upon the straight 
line decay and increase in amplitude with echo number 

In this case the interference was not due to mechanical 
misalignment which gives similar graphs but to the development of 
a crack at the edge of the transducer. As the aluminum exposed 
by the crack reacted with the water and was eaten away, the 
electrostatically driven surface was no longer circular and the 
piston like behaviour was impaired. The interference effects 
became worse with time as can be seen by comparing the 1725 and 
1753 curves beyond the tenth echo. 

c. Diffraction Effects - Comparison with Theory 

Seki, Granato and Truell [ 23 ] have used a pulse echo 
system to investigate diffraction effects in the ultrasonic field 
of a quartz crystal acting as a piston source and verify William's 
[ 27 J theory. Our pulse echo system differs in that the pulse is 
repeatedly reflected from a reflector the same size as the source. 
At first thought, the method appears equivalent (assuming perfect 
reflection) to the direct radiation from a piston. In fact, the 
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physics more closely approachs the condition of a sound beam that 
is passed through a series of apertures. At each aperture there 
will be a new diffraction (near) field. This problem has not 
been solved, to our knowledge. 

The comparison with William's theory is solely for interest. 
It must be emphasised that the determination of attenuation by 
subtracting the calibration-slope from the sea-data-slope bypasses 
the problem of diffraction that is being discussed here. 

To compare our field pattern with Seki's we have started 
with his data for ka = 50 (curve a). This curve was then modified 
to include the loss due to imperfect reflection at the transducer 
and reflector (Section I) and yielded curve b. It was assumed 
that this loss was independent of frequency. 

The comparison experimental data used was that of Section 
H, calibration of the system, since at the frequency of interest 
the theoretical absorption by fresh water was negligible. 

The active radius of the transducer was 12 cms so that, 
at 96 Kcps 



ka - 2 -tt. f . a 
c 



- 49 

This is so close to 50 that the variation in diffraction 
field due to it is negligible. 

It was also determined that at 96 Kcps a 2 /2v_was 93 • 5 cms 
so that the distance travelled by the sound waves between echoes 
was I .64 a^/^_. Curve c was then drawn with echo amplitude plotted 
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against distance travelled by sound in units of a • Different 
starting heights are used for a) and c) to preserve clarity. 

When d, the distance to the first pulse of the straight 
line section of the graph was plotted against a^/^ for the 
frequencies between 24 and 240 Kcps the straight line of Fig. 3.4 
was obtained. Its equation was 

d = 1.25 + 13.8 a2/^ 

This distance is the range beyond which the attenuation in db is 
proportional to the distance travelled. 

The dashed line for 2a also shows that diffraction 
effects and significant divergence losses occurred for beyond the 
nominal end of the near field. 

Diffraction effects varied between transducers and changes 
have been observed in a single transducer as shown by Fig. 3.5. 

The additional noise in photograph 3.5a is due to a wide band 
width setting on the electronic filter . Despite the changes in 
early pulse amplitudes, graphs of pulse amplitude against echo 
number yielded straight lines which were of the same slope within 
the limits of accuracy of the experiment (0.06 db/m). 
d. System Alignment . 

To prevent interference effects of the type shown in Fig. 

3.2 it was necessary to align the transducer and reflector so that 
their faces were as nearly parallel as possible. As shown in Fig. 
3.6 the transducer mount could be moved along the direction labelled 
"h" and rotated about its supporting arm, while the reflector was 
mounted on three adjusting bolts so that its direction might be 
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adjusted both vertically and horizontally. 

The initial mechanical alignment was made by adjusting 
the transducer mount and the reflector until the distances "s" and 
"d" were the same between all corresponding points on the 
circumferences of the transducer and reflector. 

The system was then placed in water and operated at 192 
Kcps when gross misalignment was indicated by a decay pattern much 

t 

more irregular than those illustrated in Fig. 3.5. Minor adjust- 
ments were then made at this frequency* the most sensitive within 
the normal operating range of the equipment. 

When a photograph had been taken and a plot of pulse 
amplitude against frequency found to yield an ultimate straight 
line as in curve c of Fig. 3.3* the system was regarded as being 
satisfactorily aligned. 

This procedure was difficult under laboratory conditions 
and proved more difficult while operating from the N.E.L. Tower* 
where the system was mounted as shown in Fig, 3.7. Since the 
transducer mount was easily rotated* despite locking screws* 
handling difficulties prevented the maintenance of the accurate 
allignment needed for higher frequency operation, when working from 
the U.S.N.P.G.S. Oceanographic Research Vessel. 

It is therefore recommended that the mounting should be 
redesigned so that it can be adjusted accurately in the laboratory 
and then transported safely to the observation site, 
e. Choice of Reflector Material 

To keep reflection loss as low as possible reflectors of 
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steel, styrofoam and corprene were made and used in succession. 

With a polarising voltage of 200 V. and an input drive of 20 V. 
peak at 60 Kcps applied to the transducer, the amplitudes of the 
first echo were: steel, 33 mv.j styrofoam, 31 mv.j corprene, 

25 mv. 

The steel reflector was therefore selected; it had the 
additional advantage of mechanical strength so that its reflection 
coefficient would not be altered by compression as its depth increased. 

f . Selection of Operating Voltage 

After the reflector had been selected the polarising voltage 
and input drive were varied and the first echo amplitude observed, 
the results being as shown in Fig. 3.8» On the basis of this graph 
it was decided to apply a polarising voltage of 450 V., the maximum 
which would be accurately set with the power supply in use, and 
to use an input drive of up to 50 V peak. This input drive was 
reduced at higher frequencies to maintain a better input pulse shape. 

g. Determination of Standard Deviation . 

In order to determine whether variations in the Held data 
were due to the medium or the system and to establish the confidence 
level for the calibration, the mean and standard deviation were 
determined at 60 Kcps and 192 Kcps. A series of photographs was 
taken as rapidly as possible to minimize changes in the tank water 
and equipment. The attenuations were then determined from each 
photograph. 
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60 Kcps Result . (27 samples) 



Mean X = 1.62. db/meter 
27 _ 

(X.! - X)^ = 1.58 x 1(T 2 (db/meter)^ 

j=l 

Variance S^ = 6.10 x 10~^ (db/meter)^ 

Standard Deviation S = 0.025 db/meter 
192 Kcps Result. (26 s ample s ) 

Mean X = 0.89 db/meter 
Standard Deviation S = 0.015 db/meter 

Fig. 3.9 shows the attenuation in db/m vs 
sample nunber at 60 heps. The sample points are distributed in 
much the same way as some of the field data, particularly that in 
which depth and frequency were held constant while photographs were 
taken at 2 minute time intervals. The noise was not significantly 
different in the field. The field curves have therefore been 
smoothed by up to twice the standard deviation for system noise, 
h. Calibration of System . 

Knowing the mean and standard deviation at 60 and 192 
Kcps, it was desirable to determine the number of photographs to be 
taken at each frequency to ensure that the error was less than 
1 % at a satisfactory confidence level. Since the standard deviation 
was somewhat worse at 192 Kcps than at 60 Kcps all calculations were 
performed at the higher frequency. 

When n is the number of observations 
6 is the permissible error 
s is the standard deviation 
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and is from standard normal tables 

JL 



■fti £ 



= k £ 



s 

so that for specific values of n, k^ can be determined and the 
confidence level found. 



NUMBER OF 


CONFIDENCE 


OBSERVATIONS 


LEVEL 


n 


$ 


5 


76 


6 


81 


8 


86 



Table 3.1 



From the results shown in the above table it was decided 
to take six photographs and establish a confidence level of better 
than 80$. 

In preparation for the calibration a sonar tank with a 
capacity of about 10,000 gals, was emptied, washed down and then 
refilled with tap water to which an algicide had been aded. After 
filtering for one day to remove particulate matter, the tank was 
allowed to stand for four days so that any bubbles of air produced 
by the filling might rise to the surface or be absorbed in the water , 

Six exposures were then made at the frequencies of interest 
and used to establish the standard values for attenuation due to 
the system in bubble^free fresh water. 

At all frequencies the pulse amplitudes for the first and 
last photographs were measured from pulse 2 until they became too 
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small for valid measurement or the last pulse on the photograph 
had been measured. These pulse amplitudes were then plotted 
against echo number as before to determine the straight line section 
of the graph after which measurements from the other four photo- 
graphs were confined to this section of the graph. 

From graphs of the form shown in Fig. 3.3, curve c it was 

possible to determine the pulse amplitude P n at the beginning of the 

straight line or extrapolated beyond it and a second pulse amplitude 

PrH-m * m ^jounces beyond P n but still large enough to be measured 

accurately. Then loss m bounces = 20 log £ n db. 

P m+n 

From this the loss in db/bounce was calculated and then converted 
to db/meter. knowing that the separation between the transducer and 
reflector was 76.4 cms. 

The results of this determination yielded the curve shown 
in Fig. 3.10. Apparently algae were still causing excess attenuation 
centered about 3,50 Kcps. So that excess attenuations could be 
determined from the field data the dashed line was drawn in and 
has been treated as the standard reference attenuation due to 
divergence and reflection loss, 
i. Reflection Losses. 

If the attenuation loss at 240 Kcps is used to measure the 
reflection coefficient of the system, divergence losses may be 
neglected since the linear section of the graph used to calculate 
attenuations corresponds to great distances from the source. 

Minimum measured attenuation =0.87 db/bounce 
where the path length for one bounce is 1.53 meters. 

At the frequency of 240 Kcps 
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Absorption due to fresh water 0.02 db/meter 

Absorption/bounce due to fresh water = 0.03 dh/meter 

Total Reflection Loss = 0.84 db/bounce 

10 log^Q —J—i =0.84 
1 2 



System Reflection Coefficient 




= 0.83 . 



Theoretical Reflection Coefficient of steel is given by 



since 



Steel 



, />2 c 2 ~ gl . 2 

p 2 C 2 + C 

= 0.94 



^ steel X transducer 




Reflection Coefficient of transducer = 0.88 . 

This determination has only been made for one transducer. 

If the reflection coefficients of other transducers differ from this, 
erroneous values of excess attenuation will be found when the 
standard attenuation is subtracted from the field determination of 
total attenuation. 

j . Fresh Water Algae Experiment 

The anechoic tank was drained, washed down and refilled and 
then allowed to stand for 5 days without filtering or chemical 
treatment so that the algae present could develop. Previous 
absorption experiments had shown evidence of bubbles after two days 
standing so that this waiting period also allowed bubbles to escape. 
Although considerable attenuation was measured between 120 



and 200 Kcps as shown in Fig. 3«H> no back-scattering (Section 4) 
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was detected nor did velocity measurements (Section 5) indicate 
any dispersion. 

A sample was taken with a plankton net and an analysis 
and count kindly performed by Dr. E. C. Haderlie of Monterey 
Peninsula College, the algae population being: 

Desmids (20-40 yu-) population = 20 x lCp/liter 

Ceratium Dinoflagellates (~90yltO population = 0.4/liter 

Volvox (lOOyu/) population = 0.2/liter 

This experiment showed that phytq-plankton were capable of 
producing considerable attenuation without detectable scattering or 
dispersion at selected frequencies. Whether there were minute 
bubbles associated with the photosynthesis of the plankton could 
not be determined. 

k. Absorption in San Diego Bay 

Attenuation photographs were obtained at known times and 
depths. Using the Precision Scientific Company Galvanic Bell 
Oxygen Analyser simultaneous observations of temperature and 
oxygen content of the water were obtained from a thermistor and 
oxygen electrode mounted on the equipment as shown in Fig. 3.7* 

The Winkler tests needed to calibrate the oxygen electrode were 
kindly performed by Mr. F. Crowe and Mr. J. T. Hart of the Scripps 
Field Annex. 

The photographs were analysed and graphs drawn 
to determine the total attenuation. From these results the 
attenuations due to the system were subtracted leaving the 
attenuation due to the medium, as shown in Figs. 3.12 to 3.17. 
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There was little change in attenuation over a 30 min. 
period when observations were made at a fixed frequency. 

When the frequency was fixed and the depth varied a 
consistent attenuation pattern was produced, which appeared to 
move slowly with time. Additional observations over 24 hr. periods 
are needed to confirm that a diurnal change was being observed. 

Small bubbles may be generated during the daylight hours 
causing the attenuation at frequencies above 96 Kcps to increase 
during the day. After sunset these bubbles would no longer be 
generated but would be absorbed into the water or coalesce into 
larger bubbles so reducing high frequency attenuation (See Fig. 

3.16) . Their coalescence into larger bubbles would increase 

attenuation below 96 Kcps. as shown in Figs 3»13 and 3»14 > where 
the evening attenuations are greater than the morning attenuations. 

At 152 Kcps and 192 Kcps there was a definite relationship 
between absorption and depth, with maximum absorption near the 
surface suggesting that the absorbers were organisms requiring 
light for photosynthesis. 

Attenuation did not correlate with either temperature or 
oxygen content of the water. However,. the water was isothermal 
(in general + 0.05°C or + 0.10°C) and the oxygen content varied 
little from 100$ saturation so that these quantities could have 
had little effect. 

1. Comparison of Attenuation Values with Published Data . 

In the experimental work a short sound pulse passed many 

times through a small volume of water so that the measured attenuation 

was that of a very short path length and was independent of depth, 
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surface and bottom effects. The data published in Physics of ^>ound 
in the Sea, Part 1, Transmission, measures loss over thousands of 
yards where the attenuation was effected by the surface and the 
bottom, temperature gradients and the relative depth of the 
transmitter and hydrophone in use. 

If Fig. 14, P102, graph of average transmission anomalies 
in isothermal water is examined the attenuation of 24 Kcps is 
4.5 db/kyd = 0.0049db/m and the attenuation at 60 Kcps is 
12.5 db/kyd - 0.014 db/m. 

Similar values of attenuation are given in Kinsler and 
Frey, "Fundamentals of Acoustics". 

From Fig. 3.12 from data at a depth of 10 ft. the average 
value of attenuation $t 24 Kcps was 1.1 db/m. about 220 times the 
accepted value and at 60 Kcps it was 0.35 db/m, only 25 times the 
accepted figure. 

Our data was obtained in an environment which would have 
been unsuitable for longer path transmission loss determinations 
because of its shallowness. If our near surface operation was made 
in water with many more bubbles and much more plankton than would 
be found in a deep water experiment abnormally high attenuations 
would certainly result. 

The only previous work giving attenuations of the same order 
of magnitude is that of Hansen and Barham [^9] who used tanks resonant 
at 7.9 Kcps in the laboratory and 2.7 Kcps in the sea to determine 
attenuation by measuring reverberation time. Working in "red water" 
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areas they found attenuations greater than the maximum measurable 
with their equipment, 0.405 db/m at both frequencies. Although 
the order of magnitude of their results agrees with ours it must 
be pointed out that 

a) their equipment measures absorption not scattering and absorption 

b) their equipment measured absorption for plankton, not for 
bubbles and plankton. 

Other possible explanations of our "high" value of 
absorption have been considered: 

1) Cavitation loss: This is rejected because our maximum S.P.L. 

was 70 db re 1 jib&r whereas cavitation would require S.P.L. of 
120 db. re l^ibar at the surface and greater than 120 db. re 
l^ibar as depth increased. 

2) There may have been undetected bubbles attached to our trans- 
ducer or reflector. This possibility cannot be completely 
rejected, although the self-consistency of our data suggests that 
it is an unlikely explanation. 
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4o Detection of Backscatter, 
a„ Theory,, 

A great deal of work has been done in the study of scat ter ers 
in a sound field,, The scattering cross-sections ( £<. ) of small 
solid or liquid particles is much less than the geometrical cross- 
section for 

kR = 2TTR/^ 

where R = Radius of the scatters?”. 

A = wave length of the soundo 
This represents the well known condition of Rayleigh 
scattering where the scattered intesity is proportional to the 
fourth power of the frequency. 

The scattering cross-section also depends ©n the density 
and elasticity of the particle and the medium^ the characteristic 
impedance of the scatter material ( j? C 4 ) must be different frcsa 
that of the medium for scattering to occur. 

Figure 4,1 illustrates the variation of Si with kR for the 
two extreme cases of heavy <, rigid sphere and an ideal air bubble. 

When kR is greater than 2 S both curves coincide and show 
oscillations about the mean falue of unity for ratic o.- '--R i , 

\ 
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Figure 4.1 
Comparison of 
Scattering and 
Geometrical Cross- 
Sectiont. 



It was assumed that only those bubbles at or near 
resonance are the cause of any significant backscatter at any 
particular frequency. This is so because the Q of the bubble is 
of the order of 10 and the effective scattering cross-section at 
resonance is enormously greater than the geometrical cross-section 
of the bubble 

4TTR 2 

S( WQ ) ^ _ S 2 (wo) m 4tT 

Cross-Sectional Area r2 §2 (wo) 



Where S is the bubble damping constant due to shear vecosity, 
thermal conductivity and resadiation as described by Devin [ 4 ] . 

This relative scattering cross-section is plotted in Figure 4.2 and 
varies from 918 to 178 at frequencies of 20 Kcps and 200 Kcps 
respectively, the frequency range of the apparatus. 
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Figure 4.5( a ) shows the region with system reverberant noise 
only while Figure 4.5(b) shows the increased amplitude detected 
when bubbles were introduced. When the noise is subtracted from 
the signal- plus- noise the difference is a measure of the back- 
scattered energy from the bubbles; this is plotted in Figure 4.6. 

The photos of Figure 4.5 serves also to indicate the 
difficulty of extracting bubble scatter evidence in the presence of 
reverberant scatter from the anechoic tank walls as well as the 
system supports. 

An experiment was conducted in the laboratory using a 1.5 
inch diamenter styrofoam sphere to determine the effect of scatter 
position on signal. The experiment showed that the signal level 
changed by no more than 10$ as the sphere was moved- throughout- the volume, 

V l. « - } 
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thus indicating the relative uniformity of the field. For this 
sphere and the frequency used the resulting ka of 4 indicated that 
the scattering cross-section should be approximately equal to the 
geometrical cross-section; the experiment verified this and confirmed 
the validity of the technique used. 

The scatterer position experiment is described in Appendix 
XI. 

c. Data Processing. 

1. Correction for Electronic Noise. 

The backscatter is to be evaluated from the amplitude 
H of Figure 4.7. First, however, it is necessary to use the average 
of the reflector echos A and B (Figure 4.7) as a measure of the 
amplitude which would be obtained if a "perfect" reflector were at 
the midpoint of the system. 

Hence : 

A + B 

— - — = R = the interpolated echo level for a perfect 
plane reflector. 

In a s imi lar manner the amplitude of H midway between the 
two echos can be associated with the amplitude of the scatter, 
that is: 

H-*’ S g ^ = the total scattering cross-section in the 

insonified volume, including the equivalent 
scattering cross-section of the system noise. 

Assuming incoherent signal and noise we have; 



H 2 = Eo 2 + N 2 
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(4.1) 



where h is t.r -ovr i -\ 44 electronic hoxse voltage 
(3.53 iav) aril * ’ itter voltage. 




riraire a., 7 nohos J and 4 *t jh he; s 
airce the area of the r^f^ector (, 0 j is Kn^wr, the 



proportional! v - c a* oe 



hScatt e ring J r w s c - i t c t i c n \ 
within lin .^r J uln..^- 

nefe^ctiur area j 
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°st is the total so it .er^nj cross-section of the 
transdacer/ref^octor volume to normalise thxs to a unit 

volune °st is divvied bp /; 
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(4.3) 



>s ^ y 



In this investigation S s is expressed in Cm^/m? „ 

2 0 Correction for Divergent Field. 

The fact that the transducer radiation pattern 
changes with frequency results in the insonified volume being a 
variable as the frequency is changed. 

Consider the situation as shown in Figure 4.8 




Figure 4.8. Divergent Field Situation. 

A pulse leaving the transducer backscatters as it sweeps 
over the truncated conical volume created by the divergence angle 0 
(for this experiment 3° < © < 9°). The fractional increase 

in volume over that of the original cylindrical volume becomes s. 
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V tt<Z*c/ e 

4¥ ^ vac/* Tan 9 
v 7Ta e c/ 



^fa»Q 

V (4.4) 

The reflected pulse approaching the transducer 
also diverges and causes forward scatter from bubbles ; the 
forward scatter from the cylindrical volume will be hidden in the 
main echo but the forward scatter from outside the cylinder will 
arrive later due to the longer path length and will appear between 
echos. This devergent region from the reflector is subtended by 
the same angle 6 thereby doubling the factional volume increases 

(4.5) 



AV 

7 



2d tan 6 
a 



The ratio of effective folume to cylinder volume 

becomes: 

V = 1 + f tan 0 (4.6) 

The angle 9 is taken to that of the half power 
half beam width and is given bys 

6 = sin” 1 0.0212 A ^ 0.0212A (4.7) 

Because backscatter actually comes from the larger 
volume 7* the scattering cross-section per unit volume requires a 
correction factor of; 



49 



2 



(4.8) 




and (4.3) becomes: 




The range of Divergent Field connection factors is 
given in Table 4.1 for frequencies between 20 Kcps and 60 Kcps. 
Table 4.1 Divergent Field Connection Factors (K^) 
Frequency (Kcps) 



20 


0.600 


30 


0.676 


40 


0.730 


50 


0.764 


60 


0.794 



3. Correction for Pulse Duration. 

Consider the two successive pulse as in Figure 4.9 




Figure 4.9 Pulse Correction Situation 
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During the time of the pulse ('fc' ) any backscattered 
energy will be obscured by the much greater level of the pulse, 
to correct for this a comparison is made between energy levels 
rather than power levels, assuming that the scatter power is 
spread over the full time between the pulses (T.) 

Then 4.8 becomes: 



S *=-f ( J V 2 -) 2 ^) 

In general for this experiment! 




( 4 . 10 ) 



4. Correction for Isotropic Scattering 

Because the scattering from any resonant bubbles 
will be isotropic in nature the transducer will only intercept a 
portion of the spherical scatter in an amount depending on the 
position of the bubble within the transducer/reflector volume. 

In order to obtain a correction factor for this 
effect it was necessary to obtain a value for the mean solid angle 
subtended at the transducer assuming a homogeneous distribution 
of bubbles within the volume. 

Consider a bubble at position P in Figure 4.10. 
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Figure 4.10 

Geometry of Isotropic Scattering Correction Factor 
The solid angle -fl is given by: 

■ft = 2 TT (l-cos<j)) (4.11) 

and 



CoS(j) -Jg + Cos 



From the Law of Cosines; 



Cos 2(f) = a 2 ± b 2 - c 2 
T 2ab 



where: a 2 = (R-r) 2 + X 2 

c 2 = 4R 2 

b 2 - X 2 + (R+r) 2 

hence: il - 2 TT (l - J l + Cos 2(T ) (4.12) 

expanding: 

Cos2(f) = X 2 + r 2 - R 2 

y/r2x2 - 27 . 2rX + X4 + 34200 - 370r 2 + r^ + 

. _J2 



+27.2 rX 2 + r 2 X^ + 3^0 X 2 



r (403) 



Cos 2(|» = 2B 
A 



n -2ir (1 -Ji+ -B-) 



(4.14) 



To find the mean solid angle fi it is necessary 
to perform a double integration in the two variables r and X: 

R L 

fi = 2 7f | f -0- rdr dx 

J J TT R2L 

r=0 X=0 



n = jl 

R2l 




R L 

fi rdr dX 

R=0 X=0 

Substituting the expression for fi 
L R 






X=0 r=0 

Substituting the values: R = 13.6 cm 

L = 76.4 cm 



(4.15) 
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76.4 13.6 



A- 



1124 



If 



drdx 



X=0 r=0 

76.4 13.6 



n- 



2TT- i_ 

1124 




/X 



drdx 



( 4 . 16 ) 



An approximate value was obtained for il by 
finding the mean value of fi along the axis, (for the case r=0). 
Hence: A X^ + 34200 + 370X 2 

B (X 2 R 2 ) \ 



Substituting appropriate values gave: 

fi ~ 0.5 IT 



To find the total scattering ross-section per 
unit volume it was necessary to multiply the experimentally 
obtained values by where 



— Total Spherical Solid Angle = = g 

3 n 0.57T 

5. Tank Calibration for Reverberation. 

For various data rums S s can be calculated and 
plotted as a function of frequency. In general, however, there was 
an unknown amount of interference due to the transducer support 
structure. An attempt to correct for this a number of calibration 
runs were made in the USNPGS ancohoic tank. It was found that most 
of the reverberent level measured in the anechoic tank was due to 
scattering from the tank itself. Orientation of the transducer- 
reflector had a substantial effect on background and so the curve 
of background reverberation obtained in the "calibration" over- 
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corrects. The magnitude of the correction curve was considerably 
less than the signal levels and any correction was less than 1C#>. 

6. Summary of Corrections for Backscatter Analysis. 

Combining all the correction factors (4.10) becomes: 



S 

s 




K1K2E0 

R 



•) 



finally: S 3 = — ■ ) (4.18) 

where K = K]_K 2 

Values for K over the frequency range 20-60 Kcps 
are given in Table 4.2. 

Tqble 4.2 Total Correction Factor (K). 

Frequency (Kcps) K 



20 


2.42 


30 


2.73 


40 


2.95 


50 


3.09 


60 


3.21 



It must be conceeded that the number and magnitude 
of the corrections that have been discussed above lend uncertainty 
to the final calculation of scattering cross-section per unit volume. 
The writer feels that the final correction itself may be in error 
by a factor perhaps as great as two. This error factor will 
reflect in all values of scattering cross-sections which follow, 
d. Analysis of Data 

An analyis is now carried out assuming that the scattering 
cross-sections obtained were due only to the presence of bubbleqj 
this permits calculation of the maximum bubble density and volume 
concentration for any particular datum. 



55 



corrects. The magnitude of the correction curve was considerably 
less than the signal levels and any correction was less than 10 $. 



6 , Summary of Corrections for Backscatter Analysis. 

Combining all the correction factors (4.10) becomes: 

s s . -fc (2&g&_) 

finally: S 3 = ) (4.13) 

where K = K X K 2 f~ K 3 

Values for K over the frequency range 20-60 Kcps 
are given in Table 4.2. 

Tqble 4.2 Total Correction Factor (K). 



Frequency (Kcps) K 

20 2.42 

30 2.73 

40 2.95 

50 3.09 

60 3.21 



It must be conceeded that the number and magnitude 
of the corrections that have been discussed above lend uncertainty 
to the final calculation of scattering cross-section per unit volume. 
The writer feels that the final correction itself may be in error 
by a factor perhaps as great as two. This error factor will 
reflect in all values of scattering cross-sections which follow, 
d. Analysis of Data 

An analyis is now carried out assuming that the scattering 
cross-sections obtained were due only to the presence of bubbleqj 
this permits calculation of the maximum bubble density and volume 
concentration for any particular datum. 
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From classical theory ^22 J , the total scattering cross- 



section per unit volume is given by: 

co 



U R 2 n(R) dR 



(4.19) 



Where n(R)dR is the number of bubbles per unit volume 
which have a radius between R and R+dr. 

An approximate evaluation of S s is carried out by 
substituting the small quantity. 



fL 



_ Wo 
Wo 



1 = R - 1 and d(R = dR 
Ro Ro 



and integrating: 



R 



S c % 4TTRo 3 n (R) 



-R 



0 H ■ ) "</'/'’ 

[(/ >-if\ ) ' 1 .! 1 - ' 



( 4 . 20 ) 



The integrand is simplified by dropping higher order terms 
of the small quantity R so that we need only consider: 

-R 



• f ; ( li 

-k 



dR 



1 2R 



2 <- 



tan 



= A 



( 4 . 21 ) 



ii; we c'jujoiiuo ^ ^u.l 

then A 10 tan -1 20R 
Amax = 15.7 
if we let R = 0.05 
A«s7.S3 



and 



a_ = 

Amax 15.7 



0.50 
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In other words, the bubbles within 5 % of the resonant 
size will cause 50$ of the observed scattering cross-section per 
unit volume (S g ). 

The 5 % bandwidth was used in the determination of the number 
of bubbles per unit volume with radii between r and r + dr (hereafter 
called "bubble density" = n(R)dR) and the ratio of bubble air 
volume to water valume (hereafter called "bubble volume 
concentration" = U(R)dR). 

The bubble density over any of the 5$ increments of 
frequency is: 

S s 

n(r)dR = 2 (4.22) 

e (wo) 

where ^(WO) is the resonant scattering cross-section of a 
single bubble as shown in Figure 4.3. 



6,(W0) = 



Hence 




(4.23) 



The bubble volume concentration over the same frequency 



increment is: 



U(R)dr = -A-TT Ro 3 n(R)dr 



(4.24) 



Data was taken in San Diego Bay using the USN-NEL 



oceanographic Research Tower as an equipment platform and in 
Monterey Bay using the USNPGS Hydrographic Research Vessel. 
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The variables considered in gathering the scattering data 
were: depth, frequency, time, day/night comparison and sea state. 

It was found that the system was frequency limited to a range of 
20 Kcps to 60 Kcps when using this technique; the oscillator had 
a lower limit of-"*' 20 Kcps and above 60 Kcps the signal to noise 
ration deteriorated to a point of unusability. 

Most of the data were gathered at a depth of 10 feet 
which seemed to be the optimum depth for scatter detection without 
being influenced by wave action and surface effects. 

In measuring the signal level it was found that there were 
interference or reverberation effects associated with the 
configuration of the apparatus; this is discussed in Section 4»E. 

The interference resulted in a standing wave type pattern 
being superimposed on the signal due to scattering; becuase of this 
the measuring technique used was to find and measure all the signal 
maxi ma and minima as frequency was varied; when all the points 
were plotted against frequency a, smoothing technique was used to 
approximate the signal due to scattering alone. Examples of the 
field data and related "smooth" curve are shown in Figure 4*11* 

The "smooth" curve is used as the basis for subsequent 

curves. 

e. Interference Problems. 

All the data gathered from the MEL Oceanographic Research 
Tower suffered from spurious scattered radiation. Leaks of 
interference were spaced at about 5 Kcps intervals; this indicated 
the possibility of an interfering source about 1 foot from the 
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latter observation might be explainable by phytoplankton 
photosynthesis during the day with emission of the smaller bubbles 
which coalesce after dark, decreasing the number of small bubbles 
and increasing the nuntoer of large bubbles which rise during the 
night and disappear with daylight again; there might also be a 
vertical movement of bubbles fixed to plankton during the night. 

A 24 hour observation at 40Kcps (figure 4.21) with readings 
every hour confirmed the increased scattering activity during the 
dark hours at this frequency; the scattering level rose sharply 
between 2000 and 0800 hours. 



"bubble" density curves for runs A to E are plotted in Figure 4.22 
to 4.24. In all cases the density distributions peak between 38 
Kcps and 47 Kcps; this is in general agreement with the work done 



the distribution of wind-initiated bubbles peaked at 45 Kcps. The 
comparison may even be better when one considers that the Glotov 
experiment was done in a tank with "clean" sea water whereas the 
investigation of this paper conserns bubbles in the ocean which 
may be attached to particulate matter which would tend to increase 
the bubble mass, reducing the resonant frequency and shifting the 
distribution peak to a lower frequency. 

There is approximately a 4:1 ratio between the highest and 
lowest density peaks observed. 

The corresponding "bubble" volume concentration for runs 
A to E are shown in Figures 4.25 to 4.27. 



Assuming that the scattering is due solely to bubbles the 



by Glotov, Kolobaev and Neuimim 




which it was observed that 
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A number of runs (F to L) were conducted to determine the 
variation of back-scatter with depth. Figures 4.28 and 4.29 show 
comparitive day/night data at 40 Kcps; in general there is a drop- 
off with increasing depth. The drop-off at night is less pronounced 
and the night scatter is less near the surface. At 50 Kcps however 
(Figure 4.30) the night scatter level is greater than the day level 
at all depths and in both cases the level is essentially unchanged 
with depth indicating that bubbles of this size (radius =64 microns) 
are more evenly distributed in the medium then the 40 Kcps bubbles 
(radius = 80 microns) which are more concentrated near the surface. 

The depth run in Monterey Bay (Figure 4.31) indicated that 
there was little change in the scattering level to a depth of 56 
feet at a frequency of 35.7 Kcps. However, there was no data 
between the surface and 10 feet where the level was probably higher. 
Another reason for the apparent lack of scattering drop off with 
depth is the fact that these data were taken in a condition of sea 
state 1 whereas most of the previous data were for sea state 2 
conditions in which it would be expected that there would be more 
activity near the surface, particularly at such a low frequency. 

Referring to the total scattering cross-sections per unit 
volume as plotted in Figures 4.18 to 4.20 the mean attenuation due 
to the "bubbles -only" assumption is approximately 1.0 db per meter 
in the frequency range 30 to 60 Kcps; this figure is obtained by 
using the expression [ 22 J ; 

a = 4 . 34 ^J-js s db/meter 



64 




66 



9tr Op 9£ OS 9? 





67 





68 



Otr 9£ Q£ 




69 



45 




70 



£5 




71 




72 





73 




74 




75 



30 




76 




77 




78 



•k. 




79 





20 25 30 35 40 45 

£>£PrH-pr. 



80 




81 




5* Velocity Measurements* 
a 0 Brief Theory* 

The basic properties of a liquid that control the velocity 
of sound through the liquid are compressibility^ fa anc * density <> 

The sound velocity can be expressed as 

0 = 1 / / <Kp (5“1) 

This equation can be derived from the equation of state 03 - 

The plane wave solution to the one-dimensional wave equation^ 

=k2 4£i_ , 

can be written 




t = Ae"^ ejwfc 

where k = 4- (5=2) 

It can be seen that will describe the attenuation of the wave with 
increasing x s and 0 will describe the propagation constant of the 
wave. 



If we substitute a complex Xj> 



)C = a + jb 



into equation 5-1 and then combine equation 5-1 with 5-2 we obtain 





For bubble free water 



(5-3) 




where e 0 is the sound velocity in bubble free water. 

Since we are looking for the effect of small bubbles in small 
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concentrations we will assume that the density of water with 
bubbles is not different from that without bubbles. With this 



in mind, if we have an equation for the compressibility of water 
with bubbles, we can calculate the velocity of sound through bubbly 
water. Meyer and Skudrzyk have derived such an equation for 
bubbles of a particular size [»] . 



K - 3V 

r 



iR?«i 2 



u> 



O 

72 



% 2 n . . . % 

■%— - 1 + j d — 

1 41 * W- 



where 

V = The volume concentration of bubbles in water = ratio of 
volume of air in bubbles to volume of water 
^0 = water density 

kt'p 1 = bubble resonant frequency 

| 

U) = frequency (independant variable) 

R = bubble radius 

d = bubble damping factor 
/ 

)\, = added compressibility due to bubbles in water. 

The real part of K. is 




( 5 - 4 ) 
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t 






The velocity in bubble water is then 



/( a+a' )p 




Fox Curley, and Larsen have obtained experimental data that qualita- 
tively confirm Meyer and Skudrzyk’s equation for a distribution 
of bubbles centered at 60 microns radius and ranging from 20 to 
120 microns radius. This bubble distribution constituted 0.02 
percent of the water volume, CO 

If we know the volume distribution of bubbles over a range 
of si z es we can integrate equation 5-4 over R (note that ^ = jj- (R)) 
and substitute the new a 1 total into equation 5-5 to determine 
velocity as a function of frequency. 

This was done using a CDC 1604 computer and the results 
for various distributions are shown in figure 5-1. The computer 
program is presented in Appendix Vii. The volume distributions put 
into the computer were divided such that volume concentration values 
were assigned to discrete bands of bubble sizes. The band width 
used was to microns. An alternative method of incrementing the 
bubble wolume distribution would be to assign volume concentration 
values to discrete bands of bubble resonant frequencies. The former 
method was used because it appears more often in the literature. 



[j, 22, 26] 
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5.2 Experimental Methods. 

In order to determine dispersion due to buboles in the ocean 
the velocimeter used must meet the following requirements. First, 
it must have continuous frequency coverage from 10 Kcps to 200 Kcpsj 
these frequencies are determined from the characteristics of dis- 
persion curves shown in figure 5-1 that are computed from the bubble 
density distribution observed by Glotov, Kolobaev and Neuimin. 01 
Second, it must be quickly and easily read in order to sweep the 
frequency range before ocean conditions change. 

Any dispersion experiment will be limited by temperature 
inhomogeneities in the ocean. Liebermann has measured these 
inhomogeneities at depths of 100-200 feet and found a typical 
temperature variation of ,05°C occurring over a correlation 
distance of 60 centimeters. This is equivalent to a velocity 
variation of 0.2 M/S or 0.01 percent occurring during an average 
period of 0.6 seconds assuming a current of two knots flows by the 
transducer. Similar data taken at the sea surface as a function 
of time is not available. Figure 5-1 shows that if the velocity 
error due to temperature is of order 0.2 M/S the smallest total 
volume concentration of bubbles in the Glotov distribution that 

could be detected in ocean velocity measurements would be about 

—8 

6 x 10~-. . Salinity inhomogeneities could account for a variation 
of 1.23 M/S per part ' per thousand change in salinity but whether 
such changes can be expected during an experiment is not known. 

Of the various methods for measuring sound velocity the time 
of flight method appears best suited for this project in that 
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frequency can be changed continuously without changing the physical 
dimensions of the transducer assembly and the same basic pulse 
echo system can be used in making attenuation and backscattering 



measurements as well. A variation on the time of flight technique, 
the "sing around" technique could not be used because it would 
require major changes in electronic systems Settings thereby pre- 
cluding simultaneous velocity, attenuation, and backscattering 



in a dispersive medium. Time of flight measurements have been made 
in the laboratory by Greenspan and Tschiegg with accuracies of one 
part in 30,000 or .003 percent. {j~\ 

Two different techniques for measuring time of flight have, 
been tried. Each will be described below. 

Envelope Detector and Echo Counter 

In this technique the time interval is measured from the start 
of the main bang to the start of a preselected echo. Figure 5-2 is 
a block diagram of the system. 

The advantage of this technique lies in its ease of operation 
and reading. The one megacycle counter gives an instantaneous 
reading of travel time. In the laboratory, travel time readings 
to the sixth echo were obtained. This provided a reading resolution 
of one microsecond in about 6300 microseconds or about 0.016 percent. 
However, the reading was sensitive to frequency and echo amplitude 
changes as well as velocity so that the accuracy of the velocity 
value was about 5%. A discussion of this error and a possible 
remedy may be found in Appendix VIII, 




Also the "sing around" technique is not accurate 
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Figure 5=2 

Block diagram of envelope detector and echo counter. 

Time Marker , 

This system provides an oscilloscope display similar to that 
of a radar range marker on an expanded scale "A" scope. Figure 5-3 
is a block diagram of the system. Figure 5-4 is a typical display 
provided by this system. Appendix IX presents a detailed description 
of the system. 

The time marker system appeared very attractive because of the 
precision with which measurements could be made. When measuring 
travel time to the 11th echo a reading resolution of one part in 
10,000 was possible. 
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The method of measurement was as follows. Time was measured 



to an easily identifiable point on the first echo. The time was 
then measured to the corresponding point on a later echo. The 
velocity was then obtained from 

c _ distance traveled from first echo to late r echo 
travel time from first echo to later echo 

Difficulty was encountered in identifying the proper point on 
the later echo that corresponded to the previously selected point 
on the first echo. It is believed that this difficulty may have 
been caused by phase shift upon reflection from the transducer and 
misalignment of the transducer reflector apparatus „ The error 
due to phase shift upon reflection was determined from travel time 
data taken in an anechoic tank filled with water believed to be 
non-dispersive^. Figure 5-5 is plotted from such data and used as 
calibration curves for data taken using the same transducers. 

It should be noted that the rapid apparant change of travel time 
with frequency at low frequencies will cause low frequency determ- 
ination of travel time to be inaccurate. This effect a due to trans- 
ition from stiffness to resistance control in the transducer is 
reflected as a phase shift which is a function of frequency.^ Any 
lack of parallelism causes extra cycles to appear in front of late 
echoes at high frequencies. This ambiguity in echo leading edge due 
to misalignment is a function of frequency which makes it difficult 
to follow the preselected point on late echoes as frequency is 
^"See Section 3H for a discussion of tank preparation. 

^See Appendix I for a discussion of this transition. 
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FREQueNtY (KCP5) 



If the measuring point on the later echo does not correspond 
to the measuring point on the first echo, the measured velocity 
plotted as function of frequency in a non-dispersive medium will 
be of the form 



c = 



d 

T 0 + JL_ 

f 



(5-6) 



where 

c = Measured velocity 
d = distance traveled 
To = correct travel time 

K = Error in selecting corresponding point on later echo 
(in cycles or franctions of cycles) 
f = signal frequency 

By fitting cUrves generated from equation 5-6 for various values 
of K' to the experimental curve, an approximate value for K can be 
obtained from the curve of best fit. The measuring point on the 
later echo can now be corrected by the amount K* Figure 5-6 ill- 
ustrates the above discussion. Lack of parallelism is not a source 
of error aice K is determined. Figure 5-7 shows how the signals 
received from all points on the reflector will always add to present 
a signal that represents an echo from an effective parallel reflector 
located some average distance from the .transducer that does not 
change with frequency. The error due to uncertainty in this distance 
will cause an error of 0.96 M/S per millimeter; this error does not 
change with frequency. 
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Reading error may be caused by instability in the delay and 
pulse width in the pulse generator^ instability in the oscilloscope 
synchronizing circuity instability in the trigger points in the 
counter and rise time of the trailing edge of the delayed pulse 
which was used as the time marker. The rise time error was 
eliminated by consistantly taking measurements with the desired 
point on the echo one centimeter below the top of the time marker 
pulse. Errors due to instability were determined by taking many 
readings at a few selected points and plotting distribution curves. 
The distribution curves are shown in figure IX-1. The reading 
error as determined from these curves is conservatively estimated 
as .03 percent. 

The field measurements of velocity are discussed in the follow- 
ing section. 
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6. Data fivaluati on 



A. a) Comparison of attenuation and back-scattering data. 

To compare the data obtained in Sections 3 and 4 it 
is convenient to convert both to the equivalent extinction cross- 
section using the relationships developed in Physics of Sound in 



where a is the attenuation measured in db/meter 

(r e is the extinction cross-section of a single bubble 
N(R) is the nunber of bubbles of required radius and 
S e is the total extinction cross-section. 

Then using the experimental data 



So = 0.23a m"l . 

To determine S e , given S s , the basic relationship for 




The basic attenuation relationship is 



a - + 4.34 cr e N(R) 



= + 4.34 S e , 



a single bubble is 



QL. = KR- 

0 Q S 



but KR = 8 r so that replacing the single bubble case by that for 



multiple bubbles, 

8 



Sc “ hr ^ 



where 8 is the total damping constant 



and 




damping constant 
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A. b) Experimental Results 



Fig. 6.1a shows the extinction cross-sections 
calculated from scattering cross-section and data obtained directly 
from attenuation measurement. In view of the small initial signal 
return from the scatterers and the large number of correction 
factors applied, agreement between the curves is satisfactory. 

Since no scattering is observed above 60 Kcps it is 
believed that the effect of microbubbles is negligible above that 
frequency. This implies that bubbles of radius less than 50 microns 
are rapidly dissolved while bubbles of greater than 50 micron radius, 
possibly attached to particles, are persistent. 

Attenuation at frequencies higher than 60 Kcps is 
probably due to some mechanism other than scattering and absorption 
by microbubbles. This is believed to be absorption by minute 
plankton with negligible scattering. 

Fig. 6.1b shows extinction 'cross-section as a function 
of frequency and indicates the frequency ranges over which the 
effects of micro bubbles and plankton are dominant. 

A. c) Conditions for Absorption without Scattering. 

Data both from the sea and from the fresh water algae 
experiment have shown considerable attenuation where no scattering 
was detected, so that bubbles cannot play a major part in this 
process at these frequencies. 

Meister and St. Laurent [ 19 ] have made a study of the 
absorption due to algae in the frequency range 15-27 mcps and con- 
cluded that scattering was not a mechanism for absorption. They 

also cone.! . i clod ltnl 1 no excess structural viscosity due to the 
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algae behaved in the same manner as the structural viscosity 
present in the liquid. 

Lax and Feshback £l3} show that for particles suffic- 
iently small with respect to a wavelength absorption will be 
dominant and scattering will be negligible. This occurs when 
ka< 0.1, where a is the particle radius 

and k = .2 Tf f 
c 

The maximum radius for absorption without scattering is 



a _ 0,1 c 

2 Tf f 

= 2,4 x IOVf microns 
where F is frequency in Kcps, 

This gives the maximum particle radius for absorption 
with negligible scattering at the frequencies shown in the table. 



Frequency 

Kcps 


a . max 
microns 


24 


1,000 


30 


800 


38 


630 


48 


500 


60 


400 


76 


320 


96 


250 


120 


200 


152 


160 


192 


125 



Table 6,1 



Both these studies confirm that plankton could give 
attenuation without scattering over the frequency range of the 
measurements, but neither suggests a mechanism which would make 
the absorption a function of frequency, although definite 
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attenuation peaks were obtained in the fresh water algae experiment, 

B, a) Bubble distribution from data . 

To determine the volume of air per volume of water, 

assuming the attenuation to be due to bubbles, use was again made 

of the N.D.R.C. relationships &0. 

If is given by d R 

R o 

where R Q is the radius of the bubble resonant at the frequency of 
interest, the fraction, p, of the total extinction cross-section 
which is due to bubbles between R - dr and R + dr, is 



P = 



— i_ tan ZBt 

j> 




2 



tan 



"1 2A. 



rr — 6 

Then the attenuation due to bubbles in this range is simply pa. 

Defining u(R) dR as the volume of air per unit volume 
for the bubbles under consideration 

pa = 1.50 x 10 + 3 u(R) db/meter 
u(R) = 6.67 x 10"^ pa 
u(R) dR = 2 u(r) dr 

= 1.33 x 10~3 pa dr . 

Bubble distribution and volume concentration bases on 
scattering data were illustrated in detail in section -4. f. but 
certain curves are repeated here for comparison with the distribution 
and concentrations obtained from the attenuation data. 
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Figure 6-2 shows bubble volume concentration obtained 
from scattering data in San Diego Bay while figure 6-3 shows the 
bubble volume concentrations needed to produce the attenuations 
shown in figure 3=12 „ The bubbles of radius less than 50 microns 
are not belived to exist in significant amounts. 

Figure 6=4 shows the bubble volume concentrations 
obtained from scattering data in Monterey Bay. 

C. Comparison of velocity and backscattering data . 

Velocity data was taken at the two locations as 
described in Appendix VI. The data were taken with the intention 
of comparing with velocity data at the same time utilizing the 
computer program of appendix VII. 

Unfortunately the DuMont 404 B pulse generator was not 
available when taking data at the N.E.L. tower in San biego. A ' 

GR 1217-A pulse generator was used instead. Its relatively 
unstable pulse width caused considerable horizontal jitter reducing 
the reading accuracy to about 0.07 percent. Figure 6=5 shows 
velocity data and simultaneous data derived from scattering data 
taken at the N.E.L. tower. The broken line is the bubble free 
velocity. It was chosen so that the computed dispersion curve best 
fits the data points. It was not possible to determine the bubble 
free velocity because simultaneous salinity readings were not taken 
for substitution into the empirical formula for sound velocity in 
sea water, [ll] The data for figure 6=5 were taken on 15 February 
1965 in mid afternoon and corresponds to figure 4=25 and 6-:„. 
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Attenuation data was not taken simultaneously with velocity data. 

Figure 6-6 is a dispersion curve computed from 
scattering data taken on the morning of 18 March 1965 in Monterey 
Bay. The points appearing on Figure 6-6 are velocity data points 
taken at the same time. Again, since salinity was not known the 
bubble free velocity can not be determined. Attenuation data as 
a function of frequency were not taken in Monterey Bay. 

Examination of figures 6-5 and 6-6 reveals that the 
velocity data do not disagree with the velocity dispersion curves 
obtained from scattering and attenuation data. Unfortunately the 
largest velocity shifts in these dispersion curves occur at 
frequencies below the low frequency limit of the velocity measuring 
apparatus. This could be remedied by using a transducer which has 
a flat response that extends to lower frequencies. It is the 
transition from stiffness to resistance control (about 4.0 Kcps 
in the present transducers) that produces a phase shift which 
prevents accurate measurement. 

Examination of figures 5-1* 6-5 and 6-6' reveals that 
the shape of the dispersion curves obtained from the Glotov bubble 
distribution and the various bubble distributions obtained in 
this paper are very ^similar. Also the velocity shift from bubble 
free velocity changes very little with changing frequency above 
80 Kcps. If we assume a Glotov distribution and take one velocity 
measurement at, say, 80 Kcps we can obtain the total bubble volume 
concentration for the entire distribution by substituting the 
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measured velocity 5 V into equation 6-1 

Vol. Cone. =' 9.71 (v-v 0 ) X 10" 8 (6-1) 

which is derived from figure 5-1 at 80 Kcps. This equation 
requires a knowledge of V c the bubble free velocity which p 
according to figure 5-1 s should be obtained below three Kcps. 
Unfortunately the instrumentation described in this paper is 
not operable below 20 Kcps. However^, this single measurement 
technique could be useful in studying the decay of bubbles in a 
wake with time by simply assuming the termianl velocity 
measurement is the bubble free velocity , V D . 
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7. Conclusions. 

It is postulated that the total attenuation measured at sea at 
frequencies below 60 Kcps is due to bubbles, the effect of the 
bubbles increasing as frequency is reduced. 

At all frequencies there will be attenuation with negligible 
scattering due to plankton of size ka<0.1. This effect is dominant 
above 60 Kcps. 

It is impossible to assess the effect of any suspended inorganic 
material in the sea since its acoustic effect would be the same as 
that of plankton (algae-type). Since the San ^iego sea water was 
very clear it is assumed that there was a negligible amount of 
suspended particles. 

The velocity shift measurements do not disagree with the above 
hypothesis . 
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APPENDIX I 



The Electrostatic Transducer 

A. Brief Theory 

Electrostatic transducers have been described in the literature 
by a number of individuals. 

To arrive at the expected forms of the receiving and transmitting 
sensitivities the following derivation is presented starting with 
the equation of motion for any point localized mechanical system 
with an applied periodic forcer 

md2x +. Rmdx + sx = F 3 wt (1-1) 

dt^ dt 

where x is the displacement, 
m is the mass. 

Rm is the mechanical resistance, 
s is the stiffness 
F is the force applied. 

If the displacement is simple harmonic then: 

X = Xe swt 

and (iV-l) beomces: 

-Xw2m + jXwRm + Xs = F (l-2) 

V 

Considering the transducer as a receiver a substitution can be 
made for the force: 

F = PA 

where P is the pressure 

A is the area of transducer face 
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and further, if a large polarizing voltage is applied such 
that the total charge is essentially constant, then the signal voltage 
developed (V) is approximately proportional to the displacement 
X: 



V = kX 

hence : 

“kVw^ra ■+• jkVwRm = kVs = PA 

— v/VX* c -l t I j i/V.'iU i I '7s — lti (1-3) 

k 

Microphone sensitivity is given by: 




so: 



-MvAn + j Mwian + Ms = L - constant 



(1-4) 



We have three u.xt indicate three distinct and different 

regions of operations 

(a) Low frequency or stiffness controlled regions 

Ms constant 
M ^Constant 

(b) High frequency or mass controlled region: 

Mw2m ss Constant 



Cons tan t 
~ W2 



In this region the response should fall off at a rate of 12 db 
per octave with increasing frequency. 

(c) Mid-frequency or resistance controlled regions 
MwRm <=« Constant 

M - Constant (1-5) 



In this region the hydrophone response should fall off at a 
rate of 6 db per octave. 

The E/S transducer with its large diameter to wavelength ratio 
is assumed to be operating in this region over the frequency range 
of the investigation (20 Kcps to 200 Kcps). 

Considering the same region, one way of expressing the 
transmitting sensitivity is: 



Si-f 

where I is the driving current for constant voltage V, 
but since: 

I = wCV 

where C is the capacitance of the transducer 

Sj, = £- but from (I— 5 ) Vft — ~~ ( constant') 

wCV w v y 



W 



so S if . Constant = constant 

Another way of expressing the transmitting sensitivity is: 

q = P 

v 7 

where V is the driving voltage for constant current. 

Since 

V = X_ 
wC 

S = but from (1-5) 'A P (Constant) 

^ w 



S y = w Constant) (i-7) 



I 

wC 



( 



constant 



) 
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In this case the transmitting sensitivity should increase with 
frequency by 6 db per octave. 

The calibration data obtained for 5 transducers is in agreement 
with the resistance-control assumption and the conclusions of 
equations (1-5), (1-6) and (1-7). 

Looking at equations (1-5 and (1-6) it can be seen that the 
impedance of the transducer should also be frequency dependent: 

from (IV-5) M = ^ constant 
r ^ 

w 

from (IV-6) Sj_ = £ pss -Sonstant • 

2 , Cons tant (1-8) 

w 

This is to be expected since the E/S transducer is no more 
than a capacitor. 

B. Design and Construction 

1. General 

Reference is made to Figure 1-1 which shows a cross-sectional 
view of the transducer. The fixed electrode is the aluminum plate 
while the movable electrode is the thin aluminum film on the outer 
side of the mylar; the one mil mylar constitutes the solid dielectric. 
The electrical circuit is completed by means of a coaxial cable 
connected to the two electrodes by means of the screw terminal and 
silver painted area shown on the diagram. A small air cushion exists 
between the fixed electrode and the mylar; the concentric grooves 
in the aluminum plate add to the air cushion and increase the 
sensitivity by adding compliance to the system. 
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By means of the air chamber in the back of the transducer 
pressure equalization is obtained which makes the sensitivity 
independent of depth. Such pressure equalization was described 
by Kuttruff and Wille (12 ); there will of course be a limiting 
depth for any given size of chamber but one of the transducers 
constructed for this investigation was successfully operated at a 
depth of 40 meters with no apparent change in sensitivity. 

The transducer diameter of 24 cm is equivalent to 8 wavelengths 
at 50 kcps; this was considered satisfactory and made for a con- 
venient handling size. 

Many hours were spend in developing techniques for the 
construction of a satisfactory transducer; so many subtle details 
were involved in the final technique that a step by step procedure 
is outlined in Appendix I-B.2. 

Once the mylar was bonded to the plexiglass rim a coating was 
required to protect the thin aluminum film from the action of the 
sea water. Such a coating has negligible effect on the transducer 
perfornance provided that its mass is not significant, its thickness 
is much less than a quarter wavelength, and it does not contain air 
bubbles. In practive a coating of Epon 828 of thickness 0.35 to 
Q.40 mm was consistently obtained; this more than satisfied the 
quarter wavelength criterion even at the upper limit frequency of 
200 Kcps. The response curves obtained for several transducers 
showed that the added mass of the coating did not effect the 
resistance-control assumption described earlier. 
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In determining the resonant frequency, resistance and 



stiffness of the system it was assumed: 

Rw S at 20 Kcps 
R = Rm + Rr 
Rm ^ Rr 
R-^ Rr 

where S is the stiffness 

Rm is the mechanical resistance 

Rr is the radiation resistance 

for ka^ 1 -*• Rr = ^oCTTa^ = 2.05 x 10^ 

hence: SpsRrw = 4»1 x 10^ . Newtons 

Meter 



The resonant frequency is given by: 




/ Mo + Mr 



where Mo is the measured nftss of the 
Mr is the added mass due to the 



radiator = 0.0228 kg. 
radiation reactance. 



let Wo = = 4.24 x 10 6 

Mo 

For a piston in the range Ka^> 1, Mr Mo 
hence the resonant frequency is : 
fo = 680 Kcps 

Front and rear views of a finished transducer are shown in 
Figures 1-2 and 1-3 respectively. 

The greatest difficulty encountered was in maintaining a 
water-tight seal for the facej if water penetrated the protective 
epoxy coating it would attack the aluminum film and gradually 
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reduce the effective) area of the transducer. An example of face 
deterioration due to water penetration is shown in Figure 1-4, a 
photograph of the face removed from transducer 4 after about a week 
of operation in the sea. 

Two different resins were used in constructing the transducers. 
Shell Epon 828 and Shell ^pon 1004. Operationally, there was little 
difference but the 828 was used in most cases because of its 
easier handling qualities. Both epons are described in Appendix 

I-B.3. 

2. Construction Procedure 
Problem Statements 

(1) To provide sufficient tension in the mylar face 
while the adhesive is drying, 

(2) To keep water from getting in under the rubber 
pressure-compensating diaphragm, 

(3) To seal the face so that the water cannot react 
with the aluminum coating on the mylar, 

(4) To keep the contained air from being forced out 
through the coaxial cable, 

(5) To provide a good electrical contact between the 
mylar aluminum film and the outer conductor of the coaxial cable. 

Construction Procedures 

(1) The following is a step by step procedure for the 
final stage construction of the transducer using Figure 1=1 diagram 
as a reference. 

(2) The area around the stuffing tube hole is coated 
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with silver paint and allowed to dry. 

(3) The coaxial cable must be fitted to the transducer 
by means of the stuffing tube. The cable is fed thru the stuffing 
tube, the shielding layed back and soldered to the body of the 
stuffing tube; the thread of the stuffing tube is coated with neoprene 
before insertion; the portion of the thread that protrudes after 
tightening must be left bare so that a good electrical connection 

can be made to the aluminum film on the face. 

(4) After the stuffing tube connection is tightened, 
the center conductor can be formed into a loop and connected to the 
aluminum plate by means of the provided screw. 

(5) Since the "connection cavity" provides a path for 
the air to be forced out it is filled with Epon; this effectively 
seals off the contained air. 

(6) Epon is spread on the back of the transducer on the 
area which will be directly under the diaphragm claming ring; this 
ensures a good seal when the diaphragm is in place. 

(7) A sheet of rubber is placed over the back of the 
transducer and the clamping ring placed in position; holes are 
punched thru the rubber when the pre-drilled holes of the ring are 
lined up with the holes in the transducer back; the rubber is 
stretched during this operation and held taunt while the ring is 
tightened down. The excess rubber can be trimmed off with a knife. 

(8) The ring area is coated with Epon to complete the 
sealing of the back. The transducer is left to dry 24 hours. 

(9) A piece of mylar out to size is placed on the 
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aluminum frame with the aluminum film facing up 5 an ohmeter should 
be used to determine which side is the aluminum side since it is 
difficult to determine this visually. 

( 10 ) The mylar is stretched over the frame by means of 
strips of electrical plastic tape adhered to the mylar and pulled 
down and over the edge of the frame. Care should be taken to keep 
the tape off that part of the mylar within the area of the frame hole 

(11) The frame is put aside and a film of Epoxy Glue 
applied to the outer ring of the transducer between the edge and 
the aluminum plate (or edge of trough if later model). 

(12) The frame is fitted over the transducer so that the 
mylar rests on the face; the edges are pressed to ensure good contact 
with the epoxy glue, the frame is then weighted to ensure tension 

in the mylar while the glue drys; this setup should be left overnight 

(13) The mylar can be roughly cut away to enable removal 
of the frame; the mylar is then "cut back" from the edge about one 
eighth of an inch leaving an outer ring of plexiglas. 

(14) Silver paint is not applied to the exposed threads 
of the stuffing tube and spread down over the front edge of the 
transducer and over a small area of the aluminum film. Making small 
groves in the edge at the area where the paint is to be applied 
helps in providing a good electrical contact. Once the silver paint 
has dried a continuity check should be made with an ohmmeter. 

(15) A mixture of Epon is prepared and spread over the 

face of the transducer , which has been placed on a steady a level 

surface j by means of a corrugated cardboard "spreader" about 4 "x 5 " 

(any material which will not scratch the aluminum film will suffice). 
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The Epon can be spread over the entire face quite uniformly with 
minimum thickness; care must be taken to ensure that there are no 
pinholes exposing the aluminum and that the coating is spread right 
to the edge of the transducer. 

(16) After the Epon has dried the capacitance of the 
transducer can be measured on a bridge; it should be ^ 0.005 uf. 

(17) The transducer is now ready for use. 

3. Epon Characteristics 

Epons 828 and 1004 were both used in transducer construction 
with curing Agent U; both were satisfactory but the Epon 828 
being a liquid was easier to use and became the perferred coating. 
The characteristics of the Epons are as follows? 

Epon 828: A pourable liquid at room temperature; widely 

used in wet lay-up laminating, casting and the 
pathing of electrical equipment, and in form- 
ulating adhesives for bonding wood, plastics, 
metal and glass to themselves and to each other; 
its low viscosity makes it attractive for 
contact pressure and vacuum bag laminating at 
room temperature. 

Epon 1104s Supplied as a solid resin and requires a solvent 
such as Acetone. Combines flexibility, adhesion 
and toughness with good chemical resistance 
properties. 

Curing Agent Us Provides fast cure at room temperature; has 
short pot life and low toxicity; gives strong 



lamination 
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Mixtures (by weight): 828 — 5 _ 

Agent U 1 

1004 = 1£ 

Agent U 1 

C. Characteristics 

1. Response 

With the help of the NEL Transducer Evaluation Center 
(TRANSDEC) calibration data was obtained for 5 of the 6 transducers 
taken to San Diego for field measurements. 

The equations developed in Appendix I-A indicated that the 
transmitting sensitivity with respect to driving current 9 in 
microbars of pressure at one meter per ampers input should be a 
constant in the resistance controlled region; Figure 1=5 shows 
such curves for the 5 transducers calibrated and with the exception 
of transducer 6 are within 1.5 db of the 89 db level throughout 
the 20 Kcps to 150 Kcps frequency range of the calibration. 

The transmitting sensitivity can also be expressed in terms of 
microbars per volt and as was shown in Appendix I-A the sensitivity 
should vary directly as the frequency or increase at a rate of 6 
db per octave; Figure 1-5 shows that such is the behavior for all 
6 transducers. 

The receiving response should,, from the developed theory ,, vary 
inversely as the frequency or fall at a rate of 6 db per octave. 
Figure 1-7 shows that the 6 db per octave behavior is followed 
reasonably well although the fall rates of transducers 5X and 6 is 
slightly different from that of the other 3 transducers in that 
the slopes of the response curves are slightly different. This 
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must be attributable to the fact that transducers 5X and 6 were faced 
with Epon 828 while 1, 2 and 3 were faced with Epon 1004. 

2. Radiation Characteristics 

The fact that the transducer diameter is quite large ( 8 
wavelengths at 50 Kcps) makes it highly directive. The directivity 
patterns for the 5 transducers calibrated are shown in Figures 1-12 
to 1-18; for all 5 were calibrated at 60 Kcps and in addition 
patterns for transducer 3 were obtained at 100 Kcps and 150 Kcps. 

In comparing the radiation behavior of the transducers with 
that of the ideal flat piston three things were considered; the 
effective radius, the db differences between the main lobe and the 
first side lobes, and the directivity index. 

The calculated effective piston radii ranged from 12.5 to 
15.5 cm with a mean of 13.2 cm; this was slightly greater than the 
actual radius of 12 cm and is not surprising since the face must be 
subject to some flexing action which results in a pattern normally 
expected from a larger piston. 

The directivity index ranged from 29.7 db to 31.6 db with a 
mean of 30.5 db; the theoretical value for a piston is 29.6 db. 

The db difference between the main lobe and the first side 
lobes proved to be much more inconsistent. The differences ranged 
from 10.6 db to 16.6 db with a mean of 13.8 db. The theoretical 
difference for a piston is 17.6 db. In seeking an explanation it 
would seem that the fact that two different epons were used in the 
face construction would be a factor, but the two transducers 
closest to the theoretical value (l6.6 db and 16.2 db) had different 
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faces, so some other non-obvious feature of the construction must 
be at fault; the thickness of the epon at the edges which varied 
from the mean thickness by as much at 15$ might be a significant 
factor. The variation in the epon thickness and the fact that the 
epoxy glue was sometimes squeezed in over the edge of the aluminum 
plate thereby distorting the circular shape of the radiator slightly 
accounts for most of the assymmetry in the radiation patterns. 

In use all the back radiation portions of the directivity 
patterns were probably suppressed by the rigid steel mounting plate. 
The mounting can be seen in Figure 1-8. 

3. Capacitance and Impedance Characteristics. 

The electrostatic transducer, like the condenser microphone, 
is a device that depends for its operation on the value of 
capacitance between a fixed plate and a stretched membrane which 
in this case is the aluminum film on the mylar. 

As in any capacitors 

Cc<l/d where d is the separation between plates. For 
the E/S transducer the dielectric is made up of the mylar and the 
air cust ion in series so that the capacitance C is given bys 

1-i- + 1_ 

C Cjj c A 

where Cj^ is the capacitance due to the mylar part of the 
dielectric. 

C A is the capacitance due to the air part of the 
dielectric. 
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finally 



C = ■_! & l M &A 

(d^Km + dm) 



(1-9) 



where Km is the dielectric constant for mylar = 3 

' £ 0 is the permittivity of free space = 8.8 5 x 10-12 

O p 

A is the area of the aluminum plate =4.53 x 10“' m^ 
d_& is the thickness of the air custion 

djj is the thickness of the mylar = 2.54 x 10“ - c m (l mil) 

9 

Using a measured value of capacitance at a depth of 10 feet 
it was found that the thickness of the air cushion at that depth 
was 1.88 x 10“5 m (0.741 mil). 

It was found that the value of capacitance increased with 
water depth due to compression of the air cushionj the increase 
was quite rapid to 5 feet and then leveled off to a linear increase 
of approximately 100 uuf per foot of depth. If the pressure 
compensation system were perfect the capacitance would not change 
at all with depth but such a system would require an exact combination 
of back diaphragm and mylar membrane stiffness relationships. 

Capacitance data taken at the NEL Tower to a depth of 45 feet 
is plotted in Figure 1=9. Curve data between the surface and 5 
feet is only approximate due to difficulties created by wave action. 

Shallow-water data was later gathered in the USNPGS anechoic 
tankj this data is plotted in Figure 1-10. 
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The impedance of the E/S transducer depends solely on the 
value of capacitance discussed above. A curve showing impedance 
variation over the frequency range of interest is plotted in 
Figure 1-11 using the average value of capacitance at 10 feet from 
Figure 1-9. 

D. Factors Affecting Sensitivity 

The receiving sensitivity of the E/S transducer depends on a 
number of factors which can be expressed as follows: 

, ^ 2 _ 

E Q a E 

Sensitivity — (1-10) 

dTV 

where Eq is the dc polarizing voltage 
a e is the effective radius 
E is the source pulse amplitude 

d is the spacing between the aluminum plate and the 
aluminum film. 

T is the tension in the stretched mylar, 

Zjf is the mechanical input impedance, 
e is the charging factor associated with the mylar 

The polarizing voltage Eq was maintained at 450 volts (upper 
limit of available power supply) throughout the investigation while 
the source pulse voltage (E) used was 50 volts peak; these values 
gave an approximate 10 si ratio and ensured that harmonic distortion 
would be minimized. 

The ratio Eg/d, the electric field intensity, will increase 
with depth as the air cushion is compressed as described earlier; 
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such action increases the sensitivity. 

In construction of the transducer pre-stretching the mylar 
determined the value of the tension Tj this was a compromise 
situation since the stiffer the mylar was the more the transducer 
behaved like an ideal flat piston and it is also a factor in the 
pressure compensation system; on the other hand increasing T 
decreases the sensitivity as shown in (1-10). 

Increasing the depth increases the mechanical input impedance 
which tends to decrease the sensitivity, however the increase 
of electric field intensity with depth is the dominant factor so 
that the overall effect of increased depth is an increase in 
sensitivity, disregarding the pressure-equalization system. 

It was found that there was a charging factor (e) associated 
with the mylar due to the polarizing voltage Eq. For a given 
polarity of Eq for a period of time would have a neutralizing 
effect so that if the polarity was again reversed there was an 
observable increase in sensitivity. This effect did not present 
any major problem and did not have any noticeable effect on the data 
as long as the polarity remained unchanged for the duration of a 
data taking period. 

In addition to al 1 of the above, the fact that the capacitance 
is depth sensitive as described earlier, means tha^ wave action 
could also be a factor affecting sensitivity, particularly if the 
transducer is operated close to the surface. As an example of this 
effect consider the followings 
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Wave velocity = 10 knots 
Wave period = 3 sec 

Wave amplitude = 6 ft (crest to trough) 

For a wave train equivalent to 20 msec the maximum change in 



depth for the period of 


the wavetrain works 


out to 


be — 0.1 ft 


the resulting sensitivity changes are tabulated in 

Table 1-1 


Table 1-1. 


Depth-ft 


C-uuf 


C-uuf 


M 


5 


15400 


45 


.0260 db 


10 


16400 


10 


.0054 db 


20 


17500 


10 


.0050 db 


40 


19600 


10 


.0044 db 



Since the transducer was always operated at depths of 10 feet 
or greater the above effect was not significant in this investigation. 
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1 db per radial division 

Effective Radius 12.5 cm Frequency 60 KCS 

1st Side Lobes Down 13*5 db Test Distance 8 meters 

Directivity Index 29.9 db Depth 6.1 meters 



Figure 1-12 Directivity Pattern-Transducer No. 1 
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Effective Radius 15.5 cm 
1st Side Lobes Down 10.6 db 
Directivity Index 31.8 db 



Frequency 60 KCS 
Test Distance 8 meters 
Depth 6 meters 



Figure 1-13 Directivity Pattern-Transducer No. 2 
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I db per radial division 



Effective Radius 12,3 cm 
1st Side Lobes Down 1 6.6 db 
Directivity Index 29.7 db 



Frequency 60 KCS 
Test Distance 6 meters 
Depth 6.1 meters 



Figure 1-14 Directivity Pat tern- Transducer No, 3 
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Frequency 100 KCS 
Test Distance 8 meters 
Depth 6.1 meters 

Figure 1-15 Directivity Pattern-Transducer No. 3 
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1 db per radial division 

Frequency 150 KCS 
Test Distance 8 meters 
Depth 6.1 meters 

Figure 1-16 Directivity Pattern-Transducer No. 3 
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144 




1 db per radial division 



Effective Radius 12,3 cm 
1st 3ide Lflbes Down 16,2 db 
Directivity Index 29,7 db 



Frequency 60 KCS 
Test Distance 8 meters 
Depth 6,1 meters 



Figure I-I7 Directivity Pattern-Transducer No, 5X 
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1 db per radial division 



Effective Radius 13.5 cm 
1st Side Lobes Down 12.2 db 
Directivity . Index 30 „6 db 



Frequency 60 KCS 
Test Distance 8 meters 
Depth 6,1 meters 



Figure 1-18 Directivity Pat tern- Transducer No,. 6 

■ % 



146 



APPENDIX II 



SIGNAL GATE 

This circuit utilizes a General Radio 1217-B Pulse Generator 
to gate an AC input signal. When used in conjunction with the 
General Radio 13 96 -A Tone Burst Generator, it provides an 
additional 50 db separation of on and off signals and it provides 
more flexibility in controlling pulse width. The circuit is 
shown in Figure II-l. 

The gate consists of two sections. The first section is a 
common emitter amplifier whose collector bias voltage is supplied 
by the negative terminals of the pulse generator. The width 
of the collector bias pulse determines the width of the resulting 
wave pulse output. The bias resistors are adjusted for maximum 
gain. Since high gain is the principal means for obtaining a high 
on/off signal ratio a high gain transistor was used for this 
amplifier. The forward DC current transfer ratio for the 2.N1380 
is listed as 

30 minimum 
hFE = 100 typical 

300 maximum 

The purpose of C-^ is to present a high impedance to the signal 
generator and provide DC isolation. The value of was determined 
by decreasing it until the peaks of the first cycle of a 200 Kcps 
wave pulse were the same amplitude as the subsequent peaks. 

The collector output voltage contains a negative five volt 
pedestal which is removed by the second section. 0- and clamp 



PULSE GEM. 
(-) TERMINAL 



i 
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SIGNAL GATE CIRCUIT DIAGRAM 
FIGURE H - 1 



the "plus" output of the General Radio 1217 -B pulse generator so 
that the min imum value of the pulse is zero at the cathode of D^. 

This positive going pulse is added to the wave pulse thereby 
nullifying the negative pedestal. R 5 provides a fine control of the 
pedestal height. The addition is accomplished by the resistive 
network consisting of R^, R 5 , R^, and R 7 . Figure II -2 shows 
typical wave pulse outputs at 40 and 200 Kcps. 

For optimum performance the collector bias pulse amplitude 
should be adjusted so that the gated signal is clipped symmetrically 
when the amplifier is overdriven. Then the oscillator amplitude 
should be adjusted for maximum output amplitude without clipping. 
These adjustments should be made every time the frequency is changed. 
When time does not permit following this procedure, the pulse 
generator should be set once at ten volts and the oscillator 
amplitude set as shown in Figure II-3. This shorter procedure will 
cost a few db in on/off signal ratio. 

Brief specifications are as follows; 

On/off ration = 6l db at 25 Kcps 
56 db at 60 Kcps 
45 db at 200 Kcps 

Minimum pulse width = 5 microseconds 

Maximum pulse width = 8 milliseconds 
Figure II-4 shows the variation of voltage transfer ratio with 
frequency. 
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APPENDIX III 



TRANSMIT-RECEIVE SWITCHING SECTION 

This section isolates the transmitter from the reciver using 
three networks. See figure III-l. 

The network consisting of D-^ and D2 offers a low series 
impedance to high voltage signals and a high impedance to low 
voltage signals. Therefore the echo signals do not see the low 
power amplifier impedance. R2 and C^ isolate the polarizing 
voltage supply from the AC signal voltages. The network consisting 
of D3 and presents a low impedance to ground for high voltage 
signals and a high impedance to ground for low voltage signals. The 
main bang therefore appears across R-^. Several volts of main bang 
to appear across D3 and which seriously saturates the filter, 
amplifiers, and oscilloscope display thus making the electronic 
switch necessary. The network at the output of the GR-1217-A pulse 
generator differentiates both pulse generator outputs and clips 
the resulting positive pulses. The remaining negative pulses 
trigger the electronic switch so that channel A is on during the 
main bang and channel B is on at all other times. The switching 
transients appearing at the output of the electronic switch are 
shown in Figure 2-4. They were not wide enough to saturate the 
receiver. The transmit-receive switching section provides about 
60 db of rejection. 
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TRAMSMlT - RECEIVE SWITCHING SECTION 



APPENDIX IV 



SYSTEM CHARACTERISTICS 

Overall characteristics of the pulse echo system as one 
unit are presented below. 

The usable frequency range is 16 Kcps to 200 Kcps. The upper 
frequency is limited by the band pass filter, and low electrical 
signal to noise ratio.. Power amplifier distortion at high 
frequencies was noticable but not serious. The low frequency is 
limited primarily by the divergence of the transducer beam. 

The received signal to noise ratio as a function of frequency 
is shown in Figure IV-1. Although the broken curve is a true 
indication of receiver performance, the solid curve is more an 
indication of system usefulness since the fourth echo was 
arbitrarily fixed as the first usable echo and it includes transducer 
divergence effects. The electronic switch was the primary offender 
in contributing noise to the system. The transmitted signal to 
noise ratio is greater than 90 db below 60 Kcps and greater than 
80 db above 60 Kcps. 

The maximum wave pulse width available for the system is limited 
by the tone burst generator to 128 cycles. The pulse repetition rate 
is controlled from the tone burst generator and may be varied from 
a continuous signal to one pulse in ten seconds. The minimum pulse 
width is five microseconds and is limited by the signal gate. 

The receiver output deviates from linearity no more than about 
0.25 db or 2.8 percent. This figure is approximately the order of the 
resolution of the instruments used in making the linearity check. 
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APPENDIX V 



EQUIPMENT SPECIFICATIONS 

The following is a list of specifications of the commercial 
equipment used in this project. Only those features pertaining 
to the project are included. 

Audio Oscillato r TS-382-DU 

Manufactured by TRAV-LER Radio Corporation 
Serial number: 1152 

Frequency ranges 20-200,000 CPS, continuously variable 
within four decade ranges 

Output signal amplitude: 0-10 volts, continuously variable 

within seven decade ranges 
Harmonic distortion: less than one percent 

Output impedance: 1000' ohms 

Tone Burst Generator 

general Radio Company 1396-A 
Serial Number: 397 

Signal frequency range: DC-500 Kcps 

Maximum input voltages + seven volts 
Input impedance: About 10 kilohms 

Ma ximum voltage output: + seven volts 

Total distortion: less than -60 db at 1 Kcps and 10 Kcps 

On/off ratio: -40 db with maximum signal input 

Output impedance: 600 ohms 

Gating voltage output: 12 volts, gate open 



-12 volts, gate closed 
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Pulse Generator 



General Radio Company 1217-B 

Pulse repetition frequency: 2.5 CPS to 500 Kcps 

internally generated. Externally controlled,, DC to 1 
MCPS. Output pulse is started by negative slope. 

Output pulse durations 100 nsec to 1 second in seven decade 
ranges 

Rise and fall times s Less than 20 nsec. 

Overshoot and noises less than 5 percent of amplitude. 

PRF stability: .05 percent 

Pulse duration stability: .03 percent 

Power Amplifier 

General Radio Company 1233-A 
Serial numbers 539 
Frequency Ranges 20 to 3 MC 
Output powers 15 watts 

Voltage outputs .1 microsecond rise time^ negligible overshoot. 
Input voltages .2 volt for full output 
Input impedances 100 K-ohms a 37 pf 

Distortions Less than three percent at rated output. 

Noises 70 db below 15 watts or .3 volt peak to peak. 

Power consumptions 120 watts at zero output 
Electronic Switch 

DuMont Laboratories type 330 

Serial numbers 1111 

Input impedances 2 megohms s 45 pf 
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Response: Flat to DC with 60 pf load, response down no more 

than 1 db at 10 MC. 

Hum and noise level: 2.5 MV maximum exclusive of s switching 

transients. 

Switching signal rise time: Two percent of square wave period 

External trigger switching rate: Zero to 100 Kcps 

Synchronizing Voltage: At least + 2 volts peark with rise. 

.02 V/ s. 

Power consumed: 160 watts 

Electronic Filter 

Spencer-Kennedy model 302 
Serial Number: 1483 

Frequency limits: 20 CPS to 200 Kcps in four decade ranges 

Rate of attenuation in rejection band; 18 db/octave per section 

Maximum attenuation: 90 db 

Input impedance: 2 megohms, 30 pf 

Maximum input signal: Six volts peak 

Output impedance: 300 ohms 

Minimum external load for operation according to specification: 
Ten K-ohms 

Noise Level: 80 db below one volt (using- detector bandwidth 

from zero to 40 Kcps) 

Power consumed; 120 watts 



158 



Amplifier 



Hewlett-Packard model 450-A 

Serial nunberss 228-06357 and 228-06379 

Gain; 20 db + l/8 db or 40 db + l/8 db 

Frequency responses 20 db gains within -g db 5 to 1 MC 

40 db gain: withinf^ db to 10 to 1 MG 
Maximum input voltage: 20 db gain: 1 volt 

40 db gains .1 volt 
Input impedances one megohm* 15 pf . 

Outputs Tenvolts maximum across 3000 ohms (or more) 
Distortions Less than one percent at rated output 
Internal impedance: less than 150 ohms 

Equivalent input noise levels 40 db gains about 40 v 

20 db gain: about 250 v 



Power consumed; 55 watts 
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APPENDIX VI 



TEST SITES 

a) The U, S. Navy Electronics Laboratory Oceanographic Research Tower , 
The Tower was chosen as the site for field work because it 
provided a stable platform ^29^ for the experimental equipment and 
was situated in water of known characteristics. 

The Tower is located one mile off Mission Beach where the 
bottom is gently sloping rippled sand and the mean tidal depth is 
56 feet. The dominant longshore current is from the north £ 30 ] 
so that water from Mission Bay or the flood control channel is 
unlikely to reach the Tower and the surrounding water is of almost 
constant salinity in the winter. 

The electronic equipment was placed on the Main Deck in the 
room at the right hand side of Figure VI-1 with the transducers, 
reflector, temperature and oxygen sensors mounted on a car beneath 
Figure VI-2, This car was firmly mounted on railway tracks and could 
be lowered to the sea bottom so that it provided a stable platform 
at any desired depth. 

The low self -noise of the Tower was an advantage which may have 
been somewhate offset by echoes from the structural members and 
marine growth shown in Figure VI-2. 

Full details of the construction and operation of the Tower 
will be given in a publication which is now being prepared by Dr, 

E, C, LaFond of NEE. 
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b) Monterey Bay 



Experiments were performed in Monterey Bay using the USNPGS 
Oceanographic Research Vessel as platform. Engines were stopped 
while measurements were being made so that effects due to wake 
and screw action were minimised. The vessel was usually moored 
300 yards off shore and 300 yards east of Monterey Municipal 
Wharf in 27 feet of water over a sandy bottom. 
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APPENDIX VII 



VELOCITY SHIFT COMPUTER PROGRAM 

This computer program computes sound velocity shift as a 
function of frequency in sea water containing a distribution of 
bubble sizes. The basis of the computation is the expression for 
added compressibility due to bubbles 



Q. - 



3V 



tOo 
CO 'I 






(5-4) 



For a distribution of discrete bubble sizes the compressibilities 
will add obtaining at frequency ( 0 . 



=. 

L 








2 . 



where i will be the index for each bubble size. The velocity at 
frequency, W, will then be; 



v(io) = 



+ a (w) 
Co** 



The damping factor as a function of bubble resonant frequency 
has been approximated by 

.322 



d ± = .027 



w i 1 °*- 
2.TT X 1(P J 
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(VII-1) 



Equation VII-1 fits Devin's curve within about five percent below 
200 Kcps. M • 

The bubble resonant frequency as a function of bubble radius 
at one atmosphere has been approximated by 




8TT x 10 6 

«i 



(VII-2) 



where R^ is in microns. 

In equation VII-2 the bubble is assumed to behave adiabatically, 
surface tension is ignored and the external pressure is taken as 
standard sea level pressure. 

The computer program provides a print out of volume concentration, 
formula coefficient, damping factor, and resonant frequency for 
each bubble sizej and a print out of sound velocity shift from 
1500 M/S, density times compressibility and frequency. The formula 
coefficient is 



Ci = 3V i /R i 2w 0 2 

The U. S. Naval Postgraduate School subroutine DRAW was utilized 
providing a graph of velocity shift versus frequency. [31] 

Figure VII-1 is a flow diagram for this program. Table VII-1 is 
a typical print out. 

Several volume distributions were put into the program. Each 
will be described below. 

Glotov . Kolobaev and Neuimin distribution . 

For 25<-R<69 y< 

V± = B(.2067Ri - 5.17) x 4.19 x 10“ 18 R? 
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For 71^^105 

V i = B x 1.91 x 10~ 13 R i 
For 107<R<325 / m. 

Vi = B x 2.57 x 10“7 /r^ 

where 

B = scaling factor 

Vi = volume concentration of all bubbles with radii one micron 
on either side of Ri« 

The above formulas are derived from information given by Glotov et. 
al. in their paper. However the authors indicate that R^ should 
appear in the denominator of the third formula. Ri was used 
instead because R^ does not cause Vi to converge rapidly enough 
causing a large discontinuity at the larget radius of 325 microns. 
Distributions from attenuation data . 

The folume concentrations derived from figure 3-12 were 
normalized to two micron increments. Two distributions were 
extracted from this data; the first utilizes all the data and the 
second utilizes the center portion of the distribution only. 

First curve: 

For 16 <• R < 26 /x. 

V ± = ( ,50R i - 8.0) x 10“ 10 
for 28< R< 60 ^ 

V ± = (■r.912R i i 30.1) x 10- 10 
for 62<R<84 /a. 

Vi = (.259Ri - 15.05) x 10- 10 
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for 86<Rc 110/*. 

V i = 4.3 x 10-VRi 3 
for 112 <■ R < 156 ^ 

Vi = (O.IRi - 7.8) x lO” 10 
for 158<-R < 416 

V ± = 3.05 x 10-3/Ri3 
Second Curve: 

For 58<.R<'84 / /< 

V ± = (.259Ri - 15.05) x lO -10 
for 86-< R<-458 

V ± = 4.3 x 10-yR i 3 
Distribution from Scattering Data 

The following formulas apply to Figure 6-2 and 6-4. 
San Diego 
For 65 ^R <149 A. 

Vi = (.283Ri “ 18 *3) x 10" 10 

for 151<> R^ 465/^ 

V ± = 0.0080/Ri3 

Monterey Bay 

For 65<R‘ < -135/^ 

V ± = (O.IRi - 6.5) x 10-10 
for 137 *-R< 465 /< 

V ± = 0.0017/Ri3 
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Figure ~vtt -l 
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PROGRAM SCAT A 

THIS PROGRAM COMPUTES VELOCITY SHIFT 
ACCORDING TO MCNEIL*S SCATTERING DATA 
CURVE A 

DIMENSION WR (20 I ) ,R( 201 ) ,C( 20 1 ) ,C(201 ) ,W(2C 1 ) ,A ( 201 ) , 

1SUMA (20 1) ,V(20 I ), ITITLEI 12) ,EE< 2C1) , VOL (20 I ) 

PRINT 800 

800 FORMAT ( 69H 1 I RESONANT FREQ. DAMPING FACTOR COEFFICIENT 
1V0LUME CCNC . // ) 

R ( I ) = RUBBLE R AC I US 

WR ( I ) = BUBBLE RESONANT FREGUENCY (RADIANS) 

0(1) = BUBBLE CAMPING FACTOR 
C ( I ) = FORMULA COEFFICIENT 

VOL ( I ) = VOLUME CCNCENTPATICN FCR BUBBLE SIZE I 
I = BUBBLE SIZE INCEX 

W(J) = FREQUENCY AS INDEPENCANT VARIABLE (RACIANS) 

A ( I ) = DENSITY X ACDEC COMPRESSIBILITY CUE TO BUBBLES OF SIZE I 
AA = DENSITY X COMPRESSIBILITY OF WATER WITHOUT BUBBLES 
SUM A { J ) = SUM OF ACDEC COMPRESSIBILITIES OVER ALL BUBBLES SIZES 
AT FREGUENCY W(J) 

V(J) = SOUND VELOCITY SHIFT FROM 1500 M/S AT FREQUENCY W(J) 

J = INOEPENDANT VARIABLE INDEX 
EE ( J ) = FREGUENCY IN KCFS 
R ( 1 ) =65 .0 
DO 50 1=1,200 

10 WR ( I ) = 8.0*3. 14 16*100000 C.C/R( I ) 

D ( I ) =0 . 027* ( WR ( IJ/6283.2 )** 0.322 
IF ( 1-42)21,2 1,22 

21 VOL( I ) = ( .283*R( I )- 18 . 3 ) *1 .0E-10 
GO TO 30 

22 VOL ( I ) = • 0080 /R ( I)*«3 

30 C( I)=3.0*VCL(I)/( (R( I ) * 1 . GE-6* WR ( I ) )**2) 

PRINT 801,(1, WR( I ),D(I )»C( I ) , VCL ( I ) ) 

801 FORMAT ( 15, 1P4E16.7) 

R ( I + 1 ) =R ( I ) + 2 • 0 

50 CONTINUE 

AA = 1 .0/2250000. C 
W( 1 )=6283.2 
DO 40 J= 1 , 200 
SUMA ( J )=0.0 
SUMV = 0 . 0 
DO TCI 1=1, 2CC 

A( I )=C( I )*( (WR ( I)/W( J ) ) **2- 1. )/( ( (WR( I )/W( J ) )**2-1. )**2 
1 + ( D ( I)«WR(1)/W(J) ) **2 ) 

SUMV = SUMV+VOL( I ) 

101 SUMA (J)=SUMA(J)+A( I) 

V ( J ) = ( 1 .O/SGRTF (AA+SUMA(J) ) J-15CC.0 
EE( J)=W( J)/6283.2 
W( J+1)=W( JJ+6283.2 
40 CONTINUE 

PRINT 400, (SUMV) 

400 FORMAT (30H0TOTAL VOLUME CONCENTRATION = 1PE16.7) 

12 PRINT 500 

500 FORMAT (30H1VELCCITY SHIFT CUE TC BUBBLES/) 

PRINT 501 

501 FORMAT ( 26H IN MCNEIL DISTRIBUTION//) 

PRINT 600 

600 FORM AT ( 54H FREQUENCY RHO X COMPRESSIBILITY VELOCITY SHIFT 

' 1//) 

PRINT 700,(EE( J ) , SUM A ( J ) , V ( J),J = 1,67) 

700 FORMAT (0PF8.2»6H KC 1PE1 7.7, CPF1 8.7, 5H M/S) 

PRINT 10CC 
1000 FORMATt 1 H 1 / / / ) 

PRINT 600 

PRINT 700, ( EE( J ) , SUMA ( J ),V( J) , J = 68, 1 34 ) 

PRINT 1000 
PRINT 600 

PRINT 700,(EE(J),SUMA(J),V( J ) , J = 1 35 , 200 ) 
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IT I TLE ( 1 ) =8HSCUNC 
IT[TLE(2)=8HVELCCITY 
ITITLE(3)=8H IN 
ITITLE(4)=8HWATER 
ITITLE(5)=0H WITH 
ITITLE{6)=8HBUB6LES 
ITITLE(7)=8H 
IT ITLE(8)-6HVERSUS 
I T I T L6 ( 9 ) =8H FREQ. 

I T I T LE ( 1 0 ) = 8 F 

I T I T L E ( I I ) = 8 H R.H. 

ITITLE(I2)=8h MARKS 

CALL DRAW(15C, EE, V, 0,0, 4H 

END 

END 



. ITITLE f C,C,4,C.2,2f9,8, 1 ,LAST) 



I 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 



RESCNANT FREQ. 

3.8665846E+C5 
3.751 1642E + 05 
3.6424348E+05 
3. 5398310E+05 
3.4428493E+05 
3.3510400E+C5 
3. 2640000E+05 
3.181567 1 E+C5 
3.1 028 1 48E+05 
3.0280482E+G5 
2 . 9 S68000E + C5 
2 . 8888276E+05 
2 .8239 101 E+05 
2.7618462E+C5 
2 . 70245 1 6E+C5 
2.6455579E+C5 
2.591C103E+C5 
2. 5386667E+C5 
2.4883960E+C5 
2.4400777E+C5 
2 . 39360G0E+C5 
2.3488598E+C5 
2.3057615E+C5 
2.2642162E+C5 
2.224141 6E+05 
2.1 854609E+C5 
2.1481 026E+C5 
2.1 1 2COOOE+05 
2. 077C909E+C5 
2.0433171E+C5 
2.0 106240E + C5 
1 .9789606E+C5 
1 .9482791E+05 
1 .9185344E+C5 
1 .8896842E+05 
1 .8616889E+05 
1 .8345109E+C5 
1 .8081 151 E+05 
1 .7e24681E+C5 
1 .7575385E+05 
1 .7 332966E + 05 
1 .7097143E+05 
1 .6867651E+C5 
1 .6644238E + 05 
1 .6426667E+05 
1 .6214710E+05 
1 .60081536+05 



DAMPING FACTOR 

1 *01731*678-01 
1 .0074674E-01 
9.9797C43E-02 
9.8883C59E-02 
9 .80C2 492E-C2 
9.71 53255E-02 
9 . 6323 44CE-0 2 
9.554 15C3E-02 
9. 477524 3E-02 
9.4033789E-02 
9.331 5 586E-0 2 
9 .26 19 38 3E-02 
9.1944C22E-02 
9.1 288433E-02 
9.065162CE-02 
9.003266CE-02 
8.9430 690E-02 
8.88449056-02 
8.8274568E-02 
8.7718964E-02 
8.7177 443E-C2 
8.6649389E-02 
8.61 34225E-02 
8 .5631 4C7E-02 
8.5140L27E-02 
8.46606C2E-C2 
8.4 1 92 C79E-02 
8.3733831E-02 
8 . 328 5 652E-02 
8 .2847 162E-02 
8.241 7S97E-02 
8.1 997E 1 6E-02 
8.1 586292E-02 
8.118311 7E-02 
8 .0787S99E-02 
8.04C0659E-02 
8 .002063 1 E-02 
7.9648264E-02 
7.92827 1 7E-02 
7.8923962E-02 
7.8571 779E-C2 
7.8225S6CE-02 
7.7886305E-02 
7.7552625E-02 
7.7224737E-02 
7 .6902 467E-02 
7.6585648E-02 



COEFFICIENT 

4. 511 9378E- 1 4 

3. 1 393589E- 1 3 
5. 8275239E- 1 3 
8.515689CE-13 

1. 12C3854E-12 
1. 3892019E-12 
1. 656C 184E-12 
1.9268349E-12 

2. 195651 4E-1 2 
2. 4644660E- 1 2 
2. 7332645 E- 12 

3. CC2 1G1CE-12 
3. 27C9 175E-12 
3. 5397340 E- 12 
3. 8C655C5E-12 
4. 077367 CE-12 
4.3461835E-12 
4. 615CC0CE-12 
4.8E36165E-12 

5. 1 526 330E- 1 2 
5. 4214 49 5 E- 12 
5. 69C266CE-12 
5. 9590826 E- 12 

6. 227859 IE-12 
6. 4567 1 56E- 1 2 
6. 7655321 E- 12 
7. 0343486E-1 2 
7.3C31651E-12 
7. 5719816E-12 
7. 84C758 IE- 12 
e. 1C56146E-12 
8. 378431 IE-12 
8. 6472476E- 12 
8.916064 IE-12 
9. 1 6488C7E- 1 2 
9.4536972E-12 
9.7225137E-12 
9.9S133C2E-12 
1. C26C 147E-1 1 
1.0528963E-1 1 

.1. 0797780 E-l 1 
1. 1 C66596E- 1 1 
1. 1 486C50E-1 1 
1. 1 C35669E-1 1 
1.060853CE-11 
1. 0203 153E-1 1 
9.8181 697E-1 2 
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VOLUME CCNC. / 

9.50CCCC0E-12 
6.6100CC0E-1 1 
1.227CCCCE-1C 
1 .753CCC0E-1C 
2. 3590C00E- 1 C 
2.925CCC0E-1C 
3.491CCC0E-1C 

4. C57CCC0E-1C 
4.6230CCCE- 1C 

5. 1 890CCCE-1C 
5.755CCC0E-10 
6.32 1CCC0E-1C 
6.8870CCOE-1C 
7.453CCCCE-1C 
8.019CCCCE-1C 
8.5850CC0E-1C 
9. 1510CC0E-1C 
9.717CCCCE-1C 
1 .0283CC0E-C5 
1 .0849CC0E-C5 
1.1415CCCE-C5 
1 . 1 981CCCE-09 
1 .2547CCCE-C5 
1.311 3CCCE-C5 
1 .3679CCCE-C5 
1 .4245CCCE-C9 
1 .4811CC0E-C9 
1 . 5377C0CE-C5 
1 .5943CC0E-C5 
1.65C5CCCE-C9 
1 .7075CCCE-C5 
1 .7641 CC0E-C5 
1.8207CCCE-C5 
1 .8773CCCE-C5 
1 .9339CC0E-C5 
1 .9905CCCE-C5 
2.0471 CCCE-C5 
2. 1C37CCCE-C5 
2. 1 6C3CCCE-C5 
2.21 69CC0E-C5 
2.2735CCCE-C5 
2.3301 CCO E-C5 
2.4184 170E-C5 
2.3235681 E-CS 
2.2336529E-C9 
2. 1482998E-C5 
2. 067240 6E-CS 
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48 1 .5806792E+05 

U9 ] . 56 10435E+C5 

50 1 .54 18896E+C5 

51 1 .5232000E+C5 

52 1 .5C49581C+05 

53 1 .UQ71479E+G5 

54 1 .4697544E+C5 

55 ] .4527630E+05 

56 1 .43616006+05 

57 1 .4 199322E + C5 

58 1 .4C40670E+C5 

59 1 .38855256+05 

60 1 .3733770E+05 

61 1 .3585297E + C5 

62 1 . 3440000E + 05 

63 1 .3297778C + C5 

64 1 .3158534E + 05 

65 1.3022176E+C5 

66 1 .28886 15E + 05 

67 1.2757766E+05 

68 1 .2629548E+C5 

69 1 .2503881E+C5 

70 1 .238C690E+05 

71 1 .2259902E + C5 

72 1 .2 141449E + C5 

73 1 .2025263E + 05 

74 1 . 1 91 1280E + C5 

75 1 . 1799437E + C5 

76 1 . 1 689674E + 05 

77 1 .1 581935E+05 

78 1 . 1 476 164E + C5 

79 1 . 1372308E+C5 

80 1 . 1 270314E+05 

§1 1 . 1 170133E+C5 

82 1 . 1C71710E + C5 

'83 1.0975022E+C5- 

84 1 .03800COE+C5 

'85 1 .0706609E+C5 

86 1 .0694809E+C5 

87 1 .C604557E+05 

88 1 .051 5016E+C5 

89 1 . 0428548E+C5 

90 1 .0342716E+05 

91 1 .0256286E+05 

92 1 .0175223E+C5 

93 1 .0093494E+05 

94 1 .0C13068E+C5 

95 9.93391 30E+C4 

96 9. 8560000E+C4 

97 9 . 7792996 E+04 

98 9.7C37838E+C4 

99 9.6294253E+C4 

100 9. 5561 977E+C4 

101 9.4840755E+04 

102 9.41 3C337E+C4 

103 9. 3430483E+C4 

104 9.2740959E+C4 

105 9 .206 1538E+04 

106 9.1 392000E+04 

107 9.07321 30E+04 

108 9.0081 720E+C4 

109 8.9440569E+C4 

110 8 .880848 1 E+04 

111 8.8185263E+04 

112 8.7570732E+04 

113 8 . 6964706E+04 

114 8.636701 0E+04 

115 8. 5777474E+04 



7.6274 120E-02 
7.596773GE-02 
7.S66633CE-02 
7 .536978 1 E-C2 
7 . 5077 9 4 5E-C2 
7 .4790694 E-C2 
7.45C7 9C2E-02 
7 .4 229 4 4 9 E -02 
7 .39552 19E-02 
7.3685 1C2E-C2 
7.34 1899CE-02 
7.3156780E-02 
7.2898372E-C2 
7.2643671E-C2 
7.2392585E-02 
7.21 45C24E-02 
7.1 9009C2E-02 
7.1660 1 36E-C2 
7.1 422646E-02 
7.11 e83S5E-02 
7.09571 87E-02 
7 .0729C7 1 E-02 
7.05C3S25E-02 
7.0281713E-02 
7.0062339E-02 
6 .9845 749E-02 
6.9631E82E-02 
6 .9420678E-02 
6 .92 12C79E-02 
6 .90C6C3 1 E-02 
6 .88024776-02 
6.8601 367E-02 
6.8402648E-02 
6.8206 2 72E-02 
6.8012 189E-02 
" 6.7820355E-02- 
6 . 76 30 7 22E-C2 
6 . 744 3 24 8E-02 
6.7257E89E-02 
6.7074603E-02 
6.6893252E-02 
6 .67 14C94E-02 
6.6536792E-02 
6. 636141 OE-02 
6.61 879C9E-02 
6.601 6 256E-02 
6.58464 16E-02 
6 .5678 356E-02 
6.55 12C43E-02 
6 .5347446E-02 
6.51 84 523E-02 
6.5023275E-02 
6.4863642E-C2 
6.47C56C6E-02 
6.4549 139E-02 
6.4394214E-02 
6 .4240 EC4E-02 
6.4088E84E-02 
6.3938428E-02 
6.378941 2E-02 
6.3641611 E-02 
6.34956C2E-02 
6.3350763E-02 
6.3207270E-C2 
6 .3065 1C3E-02 
6.2924240E-02 
6.2784659E-02 
„6. 264634 IE-02 



9.45221 35E- 1 2 
9 . 1C44 1 12E-12 
8.7733747E-12 
e.458 19426-12 
8. 1579319E-02 
7. 8717 153E-12 
7.5967326E-12 
7. 32822776-12 
7. 08949576-12 
6. 85 1 879 1 E- 1 2 
6. 6247638 E- 12 
6. 4C7E 76CE-12 
6. 1997792E-12 
6. 0C067 1 2E- 1 2 
5.81C3818E-12 
5. 6278703E-12 
5. 4529233E-12 
5.265153CE-12 
5. 124 1951E-12 
4 . 9697C73E-1 2 
4 . 82 1 3 678E- 1 2 
4.67EE737E-12 
4.54 19298E- 1 2 
4.41 02977E- 1 2 
4 . 283694 IE-12 
4. 16189C3E-12 
4.0446610E-12 
3.921 7933E-12 
3. 823C864E-12 
3. 71 63502E-12 
3.6 174C51E-12 
3. 52CC81 IE-12 
3. 4262 1 73 E- 12 
3. 33566 12E-12 
3. 2482685E-12 
3. 1639022E-12- 
3. 0824323 E- 12 
3. 0C37358E-12 
2. 9276955E-12 
2. 8E42C03E-12 
2. 783 1 447E- 1 2 
2. 7144282E-12 
2. 6479554E-12 
2. 5836355E-12 
2. 5213819E-12 
2.46 1 1 124E-12 
2. 4C27486E-12 
2. 3462 157E- 12 
2. 29 1 4424E- 1 2 
2. 2383610E-12 
2. 1 869C64E-12 
2. 1370 17CE-12 
2. 0886335 E- 12 
2. 04 16997E-12 
1.9961616E-12 
1.9519679E-12 
1. 9C9C690E-12 
1. 8674 1806— 1 2 
1. 8269700E-12 
1. 7876817 E- 12 
1-74951 18E— 12 
1.7 124209E-12 
•1 • 676 37 1 2E- 1 2 
1. 64132626-12 
1 . 6C725 1 4E-1 2 
1.5741132E-12 
1. 541 8798E-12 
1. 51C52C5E-12 



1 .9902C87E-C9 
1.91 69569E-C9 
1 .84725646-09 
1 .78089436-09 
1.71 76733E-C9 
1 .6574C97E-09 
1 .5999325E-C9 
1 .545C825E-C9 
1 • 4927 1 1 46-C9 
1 .4426EC6E-C9 
1 .39486C9E-C9 
1 .3491 31 4E-C9 
1 .3053793E-C9 
1 .26349876-09 
1 .22339C7E-C9 
1 . 1 849624E-C9 
1 . 1481269E-C9 
1 . 11 28C24 E-C9 
1 .0789 123E-C9 
1 • 0463845E-C9 
1.0151513 E-C9 
9.8514876E-1C 
9.5631699E-1C 
9.285994 1 E-1 C 
9. 0194 27 IE- 1C 
8.7629660E-10 
8.5161366E-1C 
8.27849C9E-1C 
8.0496057 E-10 
7.829C81C6-1C 
7.616538SE-1C 
7.41 162C4E-10 
7.2139877E-1C 
7.0233 1966-10 
6.8393 120E-1C 
6 . 66 1 6 7 6 5 E— 1 0 
6.4901 397E-1C 
6.324442 1 E-1C 
6. 1643374E-1C 
6.0095915E-1C 
5. 8599820 E- 1 C 
5.7152977E-1C 
5. 5753375 E- 1C 
5. 4399 1 C2E-1 0 
5.3088338E-1C 
5. 1819348E- 10 
5.0590483E-1C 
4.9400 168E-1C 
4 . 8246 9C4E-1 0 
4.7129260E-1C 
4. 6045871 E- 1C 
4.4995436E-1C 
4.397671 IE-10 
4.29885C7E-1C 
4 . 2029 692E- 1 C 
4.1099 179E-1C 
4. 0195932 E- 1C 
3.9318962E-1C 
3.8467218E-1C 
3. 7640C92E- 10 
3. 68364 1 7 6 — 1 C 
3.6055459E-1C 
3. 5296422E- 1 C 
3.4558542E-1C 
3.384 1 C86E-1 C 
3.3143354E-1C 
3.2464672E-1C 
3. 1 8C4394E-1 0 
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1 1 6 8.5195932E+04 

117 8.4 £ 222225 + 04 

118 8 .4056 1 87E+C4 

119 8 . 5497674E+C4 

120 8 . 2946535E+C4 

121 8 .2402623E+C4 

122 8. 1865798G+04 

123 8.1 335922E+04 

124 8. 061286 2E+C4 

125 8.C296486E+C4 

126 7 . 978666 7E+C4 

127 7 . 928328 1 E+G4 

128 7.87862C7E+C4 

129 7.8295327E+04 

130 7.7810526E+C4 

131 7.7331692E+C4 

132 7.6858716E+04 

133 7.65914896+04 

134 7. 5929909E+C4 

135 7.5473874E+C4 

136 7 .50232846+04 

137 7 .45780426+04 

138 7.4 1 380536+04 

139 7.37032266+04 

140 7.32734696+04 

141 7.28406966+04 

142 7 . 242881 8E+C4 

143 7 . 20 1 37546+04 

144 7. 16C3419E+04 

145 7.11977346+04 

146 7 . 07966206+04 

147 7 . 040C00CE+C4 

148 7 . CC07799E+04 

.149 6.96199456+04 

150 6.92363646+04 

151 6.88569866+04 

152 6.8481744E+04 

153 6.8 1 10569E+04 

154 6.7743396E+C4 

155 6.738C16 16+04 

156 6 . 7C2080CE+04 

157 6 . 6665252E+04 

158 6.631 34566+04 

159 6.5965354E+04 

160 6 . 562C888E+04 

161 6 . 52800006+04 

162 6.49426366+04 

163 6.4608740E+04 

164 6.4278261 E+04 

165 6.59511456+04 

166 6.36273426+04 

167 6.330680 1E+04 

168 6 . 2989474E+04 

169 6.26753126+04 

170 6.23642686+04 

-171 -6-.2056296E+C4 

172 6.1751351E+04 

173 6. 1 449389E+C4 

174 6 . 1 15C365E+C4 

175 ' 6 . 0854237E+04 



6 .25092656-02 
6.2373412E-02 
6 .22387646-02 
6.2 1052C 16-02 
6 . 1 973CC5E-02 
6.184] £596-02 
6.1711 £446-02 
6. 1582945E-02 
6.14551436-02 
6 . 1328424E-02 
6.12027716-02 
6. 10781686-02 
6 .09545996-02 
6.0832C51E-02 
6 .07 105C8E-C2 
6.05899566-02 
6.047038 IE-02 
6 .03517696-02 
6.0224 1C6E-02 
6.01172796-02 
6 .000 1 57 56-02 
5.9886682E-02 
5.97726876-02 
5.96595776-02 
5.95472416-02 
5.9435967E-02 
5.93254436-02 
5.9215 75eE-02 
5.91069016-02 
5 .8998 8 6 1 E-C 2 
5.88916276-02 
5.8785 ie9E-02 
5.86795366-02 
5 .85746586-02 
5.8470546E-02 
5.83671896-02 
5.82645786-02 
5.81 627046-02 
5.80615586-02 
5.7961 1296-02 
5.78614 1 0E-02 
5.77623926-02 
5.7664C65E-02 
5.75664226-02 
5.74694556-02 
5.73731546-02 
5.72775136-02 
5.71825236-02 
5.70881766-02 
5.6994465E-02 
5.6901 382E-02 
5.68C8921E-02 
5. 67170736-02 
'5.66258326-02 
5.6535 190E-02 
--5.6445 14 IE-02 
5.63556786-02 
5 . 6266795E-02 
5.61 78484E-02 
5 .60907396-02 



1.48CCC59Erl2 
1. 4503077E-12 
1. 42139876-12 
1. 393253 1 E- 1 2 
1 . 3656456E-1 2 
1.339 1 5236- 1 2 
1. 31315016—12 
1. 2878 167E-12 
1. 263 13C8E-12 
1. 239C718E-12 
1. 21562006-12 
1. 1 9275636-12 
1. 17C4623E-12 
1. 1487206E-12 
1. 12751406-12 
1. 1C68261 E-12 
1.08664146-12 
1 . 0669444E-1 2 
1. 0477207E-12 
1.02895606-12 
1.0 1C6367E-12 
9. 92749766-13 
9. 7528242E-13 
9.56222486-13 
9. 4 1558126-13 
9. 2 527794E- 1 3 
9. 0937C93E-13 
8. 93826476-13 
e. 78634286-13 
8. 63784446-13 . 
8. 49267366-13 
8. 35C73776-1 3 
6.21194716-13 
8.0762 1 53E- 1 3 
7. 9434582E-13 
7. 8 1 35950E-13 
7. 68654726-13 
7. 5622389E-13 
7. 44059676-13 
7.32154946-13 
7. 2C5C283E-1 3 
7. 09C9666E-13 
6. 9792999E- 13 
6. 8699656E-13 
6. 7 629C32E-1 3 
6. 6580538E-13 
6. 5553607E- 1 3 
6.4547687E-13 
6. 3562244E- 1 3 
6. 2596758E-13 
6. 1 650727 E- 1 3 
6. 0723664E-1 3 
5. 98 1 5C96E-1 3 
5. 8924563E- 1 3 
5. 8C51621E-13 
5. 7 155837E- 1 3 
5. 6356792E- 1 3 
5. 5534078E- 13 
5.4727300E-13 
5.3936075E-13 



3. 1161 9C IE- 10 
3. 0536597E-1 C 
2. 9927912 E- 1C 
2. 9335298 E- 1C 
2.8758227E-1C 
2.8196193 E- 1C 
2.76487C9E-1C 
2.71 153C8E-1C 
2.659554CE-1C 
2. 6088972E-1 C 
2.5595 188E-1C 
2.51 13787E-1C 
2.4644382E-1C 
2. 4186 604 E- 1C 
2.374CC93E-1C 
2 . 3304 5C6E- 1 C 
2.287951 1E-1C 
2 . 2464766E- 1 C 
2.206CC25E-10 
2. 1664930E-1C 
2. 1279213E-1C 
2. 09025996- 1 C 
2.053482CE-1C 
2.0175618E-1C 
1 .9824746E-1C 
1.9481 962E-1C 
1 .9147C36E-10 
1.88197446-10 
1 .84998686-10 
1 .81872C1E-1C 
1 .78815406-10 
1 .75826916-10 
1.72904646- 1C 
1 . 7004 6776- 1 C 
1 .67251536-10 
1.6451 723E- 1 C 
1 .6184221E-1C 
1 . 59224866- 1 C 
1 .56663666-10 
1.541 57C9E-1 C 
1 • 5 1 7C37CE- 1 C 
1 . 49302 1 1 6- 1 C 
1 .4695C94E-10 
1.4464e87E-lC 
1 . 42394656-1 C 
1 .40187026-10 
1 .3802479 E- 1 0 
1 . 3590680E- 1 C 
1.33831926-10 
1 .31799076-10 
1 .29807186-10 
1 .2785522E-1C 
1 .259422 1 E- 1 C 
1 . 2406 7 1 7 E- 1 0 
1 .22229176-10 
V.-2C42729E-1C- 
1 . 1866C66E-1C 
1.16928426-10 
1 . 1 522972E-10 
1. 1356370E-1C 
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] 76 6 . 0560964E + 04 

177 6.C270504E+C4 

5 . 99828 1 6E + C4 
179 5 . 9697862E + 04 

§0 |.9415603E+04 

181 5.91 36000E + C4 

182 5. 8G59016E+C4 

183 5.858461 5E+04 

184 5. 8312761 E + C4 

|5 5.8C43418E+C4 

186 5.7776552E+C4 

]87 5.7512128E+C4 

188 5.72501 14E+C4 

189 5.6990476E+C4 

190 5.67331 83E+C4 

100 5.6478202E+04 

101 5.6225503E+C4 

193 5. 5975056E+C4 

194 5.5726829E+04 

195 5.5480795E+C4 

196 5 . 5236923E+C4 

197 5. 4995 1 86E+C4 

198 5 . 4755556E + 04 

199 5.451 8004E+04 

200 5.4282505E+04 



5 .6003 5 55E-02 
5.59 1 6924E-02 
5.5830E4 IE-02 
5.S745299E-02 
5 .5660292E-02 
5.55756 16E-02 
5.5491663E-02 
5.54C8429E-02 
5.53255C7E-02 
5 . 5243C92E-02 
5.5161 175E-02 
5 .50797 62E-C2 
5 .4998E37E-02 
5.491 8397E-02 
5 .4838439E-02 
5.4758956E-02 
5.4679944E-02 
5.4601 398E-02 
5.4523 3 1 3E-02 
5.4445684E-02 
5.4368508E-02 
5.4291778E-C2 
5.42 1549 IE-02 
5.4139 d 42E-02 
5.40642286-02 



TOTAL VOLUME CONCENTRATION = 1.3114546E-C7 



5. 3 1 6CC28E- 1 3 
5. 2398798E- 13 
5. 1652C33E-13 
5. C91939lErl3 
5. 02CC540E-13 
4. 9495 156E-13 
4. 88C2927E-13 
4. 8 12354SE-13 
4. 74567 16E- 13 
4. 68C2 15CE-13 
4. 6 159566E-13 
4. 5526693E-13 
4. 49C9263E- 13 
4. 43C 1C2CE-13 
*4.370371 IE-13 

4. 3 1 1 7092 E- 1 3 
4. 254C924E-13 
4. 1 974977E-1 3 
4. 14 19C24E-13 
4. 0872847E-13 
4.0336 230E— 13 
3. 98C8965E-1 3 
3. 929C851E-13 
3.8781689E-13 
3. 82E 1286E-13 



1. 1 192979E-1C 

1 • 042133 1 E-l 0 
1.0275580E-1C 
1.0132535E-10 
9.9921323E-1 1 
9.8543 1 18E-1 1 
9.719C 142E-1 1 
9. 5861 82 1 E-ll 
9.4557597 E-l 1 
9.3276925E-1 1 
9.2019275E-1 1 
9.0784 134E-1 1 
8.957C999E-1 1 
8.8379383E-1 1 
8.72C8810E-1 1 
8. 6058819E- 1 1 
8.4928958E-1 1 
8.3818790E-1 1 
8.2727887E-11 
8.1 655E32E- 1 1 
8. 06C2222E- 1 1 
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VELOCITY SHIFT CUE TC GU8BLES 
IN MCNEIL DISTRIBUTION 



FREQUENCY RHO X COMPRESSIBILITY VELOCITY SHIFT 

/ 



1.00 


KC 


1 .58C5603E-1 2 


-. 0C26672 


M/S 


2.00 


KC 


6.4 1 393 58 E- 1 2 


-.0108234 


M/S 


3.00 


KC 


1 .4796904 E-l 1 


-. 0249692 


M/S 


4.00 


KC 


2.7 3C9566E- 1 1 


-.0460828 


M/S 


5.00 


KC 


4.5CC6175E-1 1 


-.0759422 


M/S 


6.00 


KC 


6.99C7125E-1 1 


-. 1179544 


M/S 


7.00 


KC 


1 .0658786E-1 C 


-. 1798347 


M/S 


8.00 


KC 


1 .7C38791 E-l C 


-.2874469 


M/S 


9.00 


KC 


2.2668272E-1C 


-.3823808 


M/S 


10.00 


KC 


2. 0659377 E- 1C 


-.3485055 


M/S 


11.00 


KC 


2 .05676 19E-1 C 


-. 3469582 


M/S 


12.00 


KC 


2. 0313350 E- 1C 


-.351 1020 


M/S 


13.00 


KC 


2.1C27664E-1C 


-.3547160 


M/S 


1L.00 


KC 


2.1 C 49 692 E- 1 C 


-. 3550874 


M/S 


15.00 


KC 


2.0777052E-1C 


-. 3504899 


M/S 


16.00 


KC 


2.0 126459E-1 C 


-.3395187 


M/S 


17.00 


KC 


1 .9C188C2E-1C 


-.3208393 


M/S 


13.00 


KC 


1 .7371231 E-1C 


-.2930536 


M/S 


19.00 


KC 


1 .5CS1119E-1C 


-.2545978 


M/S 


20.00 


KC 


1 .2069606E-1C 


-.2036332 


M/S 


21.00 


KC 


8. 17 5004 1 E- 1 1 


-. 1 37 9004 


M/S 


22.00 


KC 


3.22S7022E-1 1 


-.0544983 


M/S 


23.00 


KC 


-2.99 15002E-1 1 


, .0504841 


M/S 


24.00 


KC 


-1 .08192C6E-1 C 


. 1 826C74 


M/S 


25.00 


KC 


-2. 0742440 E- 1C 


.3501512 


M/S 


26.00 


KC 


-3.34311 18E-1C 


'. 5644686 


M/S 


27.00 


KC 


-4.86241 15E-1C 


.8212C58 


M/S 


28.00 


KC 


-6 . 20 75042E- 1 C 


1.0486149 


M/S 


29.00 


KC 


-7.18 528525- 1 C 


1 . 2 1 3 9 8 9 C 


M/S 


30.00 


KC 


-7.89 12450E-1C 


1.3334235 


M/S 


31.00 


KC 


-8.4C72293E-1C 


1.42G7359 


M/S 


32.00 


KC 


-8.7863755E-1C 


1 .4849029 


M/S 


33.00 


KC 


-9. 0640521 E- 1C 


1.5319023 


M/S 


• 34.00 


KC 


-9.2647964E-1 0 


1. 5658830 


M/S 


35.00 


KC 


-9.4C62358E-1 C 


1. 5898263 


M/S 


36.00 


KC 


-9 . 50 1 3636E- 1 C 


1.6059304 


M/S 


37.00 


KC 


-9.56CC018E-1C 


1.6 158575 


M/S 


38.00 


KC 


-9.5896551 E-l C 


1.6208778 


M/S 


39.00 


KC 


-9.5961931 E- 1C 


1.6219846 


M/S 


40.00 


KC 


-9.5841798E-1C 


1.61995C8 


M/S 


41.00 


KC 


-9 . 5 5 72 699 E-1 C 


1. 6153950 


M/S 


42.00 


KC 


-9 .51838 48 E- 1C 


1.6088120 


M/S 


43.00 


KC 


-9.4690572E-1C 


1.6005967 


M/S 


44.00 


KC 


-9 .41 35676E- 1 C 


1 • 59 1 C 67 4 


M/S 


45.00 


KC 


-9.351C4G2E-1C 


1. 58C4825 


M/S 


46.00 


KC 


-9.28351 1 1 E-l C 


1. 5690510 


M/S 


47.00 


KC 


-9.21 198 15E-1 C 


1.5569426 


M/S 


48.00 


KC 


-9 . 1 3 72639E- 1 C 


1.5442949 


M/S 


49.00 


KC 


-9.06CC184E-1C 


1.5312196 


M/S 


50.00 


KC 


-8.9807696E-1C 


1. 5178C55 


M/S 


51.00 


KC 


-8 . 8999037E- 1 C 


1.5041 181 


M/S 


52.00 


KC 


-8 . 8 1 76999E- 1 C 


1.49C2C46 


M/S 


53.00 


KC 


-8. 7344152 E- 1C 


1.4761086 


M/S 


54.00 


KC 


-8.65C2427E-1C 


1. 4618627 


M/S 


55.00 


KC 


-8.5651873E-1C 


1.4474678 


M/S 


56.00 


KC 


-8. 479277 IE- 1C 


1.4329287 


M/S 


57.00 


KC 


-8.3925645E-1C 


1.4182542 


M/S 


58.00 


KC 


-8. 3048841 E- 1C 


1.4C34163 


M/S 


59.00 


KC 


-8.21631 18E-1C 


1.3884280 


M/S 


60.00 


KC 


-8 • 1 269672E-1 C 


1. 3733094 


M/S 


61.00 


KC 


-8 .03752 12E-1 C 


1.3581741 


M/S 


62.00 


KC 


-7.9495245E-1C 


1. 3432845 


M/S 


63.00 


KC 


-7.8652497E-1 C 


1.329C251 


M/S 


64.00 


KC 


-7 .78 6 1 964E- 1 C 


1. 3 156495 


M/S 


65.00 


KC 


-7.71 29255E-1 C 


1.3032527 


M/S 


66.00 


KC 


-7.6453344E-1C 


1.2918171 


M/S 


67.00 


KC 


-7. 5829973 E- 1C 


1.2812706 


M/S 
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FREQUENCY RHO X COMPRESS IB IL ITY VELOCITY SHIFT 



68.00 


KC 


-7.5253968E-1 C 


1.2715257 


M/S 


69.00 


KC 


-7.472C296E-1C 


1. 2624971 


y/s 


70.00 


KC 


-7.4224420E-1 C 


1.2541CE1 


y/s 


71.00 


KC 


-7 . 37 62367E- 1 C 


1.2462915 


y/s 


72.00 


KC 


-7.333C692E-1C 


1.2289888 


y/s 


73.00 


KC 


-7 . 29264 1 1 E-l G 


1.2221497 


y/s 


74.00 


KC 


-7.2546934F.-1C 


1.2257303 


M/S 


75.00 


KC 


-7 . 2 1 9C0C0E- 1 C 


1.2196922 


y/s 


76.00 


KC 


-7 . 1 853630E- 1 C 


1.2140022 


y/s 


77.00 


KC 


-7 . 1 5 36079 E- 1 C 


1. 2086205 


M/S 


78.00 


KC 


-7.12 35805 E- 1 C 


1.2C35512 


M/S 


79.00 


KC 


-7. 0951439 E- 1C 


1. 1987410 


M/S 


80.00 


KC 


-7 .068 1 755E-1 C 


1. 1941791 


M/S 


81.00 


KC 


-7.04 25657E-1 C 


1. 1 898472 


M/S 


82.00 


KC 


-7.C182160E-1C 


1.1 857284 


y/s 


83.00 


KC 


-6.9950373E-1C 


1. 1818C77 


y/s 


84.00 


KC 


-6.9729491 E-l C 


1. 178C716 


y/s 


85.00 


KC 


-6.95 18782E-1 C 


1. 1745C74 


M/S 


86.00 


KC 


-6.93 17579E- 1 C 


1. 171 1C42 


M/S 


87.00 


KC 


-6.91 25275E- 1 C 


1. 1678515 


y/s 


88.00 


KC 


-6 .894131 3E-1 C 


1. 1647299 


y/s 


89.00 


KC 


-6 .87651 82E-1 C 


1. 16176C7 


y/s 


90.00 


KC 


-6. 8596414 E- 1C 


1. 1589C61 


M/S 


91.00 


KC 


-6 . 8434574 E- 1C 


1. 1561688 


y/s 


92.00 


KC 


-6.8279266E-1C 


1. 1 535419 


M/S 


93.00 


KC 


-6 .8 1 3C1 20E-1 C 


1. 151C192 


M/S 


94.00 


KC 


-6 . 7986794 E- 1C 


1. 1485951 


y/s 


95.00 


KC 


-6. 7848972 E-l C 


1. 1462640 


y/s 


96.00 


KC 


-6 .77 16359 E- 1 C 


1. 144C210 


M/S 


97.00 


KC 


-6 .7588 683 E- 1 C 


1. 1418615 


M/S 


98.00 


KC 


-6.7465687E-1C 


1. 12978 12 


M/S 


99.00 


KC 


-6.7347134E-1C 


1. 1 277761 


y/s 


100.00 


KC 


-6.7232802E-1C 


1. 1358423 


M/S 


.101.00 


KC 


-6.7 1 22482E-1 C 


1. 1229765 


M/S 


102.00 


KC 


-6.7C15982E-1C 


1. 122 1752 


y/s 


103.00 


KC 


-6.69121 17E-1C 


1. 1 204 354 


y/s 


104.00 


KC 


-6.68 13717E-1C 


1. 1287543 


M/S 


105.00 


KC 


-6. 67176 22 E- 1C 


1.1271290 


M/S 


106.00 


KC 


-6 . 6624680E- 1 C 


1. 1255571 


M/S 


107.00 


KC 


’ -6.6534749E-1C 


1. 1 24C261 


y/s 


108.00 


KC 


-6.6447 695 E-l C 


1. 1225627 


y/s 


109.00 


KC 


-6.6363392E-1C 


1.1211379 


M/S 


110.00 


KC 


-6.6281720E-1 C 


1. 1 197566 


M/S 


111.00 


KC 


-6 . 62C2568E-1 C 


1. 1 184179 


y/s 


112.00 


KC 


-6 . 6 1 25829 E- 1 C 


1. 11712C1 


M/S 


113.00 


KC 


-6. 6051402 E- 1C 


1. 1158613 


M/S 


114.00 


KC 


-6 . 5979 1 92E- 1 C 


1. 1 146400 


M/S 


115.00 


KC 


-6.59C9 1 C9E-1 C 


1. 1 124548 


y/s 


116.00 


KC 


-6 . 584 1 067E- 1 C 


1. 1 123040 


M/S 


117.00 


KC 


-6.5774985E-1C 


1.11 11864 


M/S 


118.00 


KC 


-6.57 10787E-1C 


1. 1 101007 


y/s 


119.00 


KC 


-6.5648399E-1C 


1. 1C90455 


y/s 


120.00 


KC 


-6 .55877 5 1 E-l C 


1. 1080198 


M/S 


121.00 


KC 


-6 . 5528779E-1 C 


1. 1C70224 


M/S 


122.00 


KC 


-6. 5471417 E- 1C 


1.1060523 


y/s 


123.00 


KC 


-6.54 15608E-1 C 


1. 1051084 


M/S 


124.00 


KC 


-6 . 5 36 1 293E-1 C 


1. 1041898 


M/S 


125.00 


KC 


-6 .53084 1 8E- 1 C 


1. 1022956 


K/S 


126.00 


KC 


-6. 5256932 E- 1C 


1. 1024249 


M/S 


127.00 


KC 


-6.5206785E-1C 


• 1. 101 5767 


M/S 


128.00 


KC 


-6 . 5 1 57929E-1 C 


1. 1C07505 


M/S 


129.00 


KC 


-6.51 10321E-1C 


1.0999453 


M/S 


130.00 


KC 


-6 . 50 639 1 5 E-l C 


1.0991605 


M/S 


131.00 


KC 


-6 .50 18672E-1 0 


1.0983954 


M/S 


132.00 


KC 


-6.4974552E-1C 


1.0976492 


M/S 


133.00 


KC 


-6.4931517E-1C 


1.0969214 


M/S 


134.00 


KC 


-6.4889531 E-l 0 


1.0962113 


M/S 
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Y 



FREQUENCY 



RHQ X COMPRESSIBILITY VELOCITY S H| I F T 



135.00 

136.00 

137.00 

138.00 

139.00 

140.00 

141.00 

142.00 

143.00 

144.00 

145.00 

146.00 

147.00 

148.00 

149.00 

150.00 

151.00 

152.00 

153. 00 

154.00 

155.00 

156.00 

157.00 

158.00 

159.00 

160.00 
161.00 
162.00 

163.00 

164.00 

165.00 

166.00 

167.00 

168.00 
•169.00 

170.00 

171.00 

172.00 

173.00 

174.00 

175.00 

176.00 

177.00 

178.00 

179.00 

180.00 
181.00 
182.00 

183.00 

184.00 

185.00 

186.00 

187.00 

188.00 

189.00 

190.00 

191.00 

192.00 

193.00 

194.00 

195.00 

196.00 

197.00 

198.00 

199.00 

200.00 



KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 



graph titled 



-6 . 48 405 59 E- 1 C 
-6.48C8569E-1C 
-6.4769523E-1C 
-6.4731406E-1C 
-6. 4694173 E- 1C 
-6.4657803E-JC 
-6. 4622267E-1 C 
-6.4587540E-1 C 
-6 .45 53597E- 1 C 
-6.45204 14E-1C 
-6.4487968E-1C 
-6. 4456238 E- 1C 
-6.4425201 E-1C 
-6. 4394838 E- 10 
-6.4365128E-1C 
-6. 43360 53 E- 1C 
-6 . 4 307 594 E- 1 C 
-6. 4279734 E- 1C 
-6.4252456E-1 C 
-6.4225744E-1C 
-6 . 4 1 9958 1 E- 1 C 
-6.41 73953E-1 C 
-6.4 1 40844E- 1 C 
-6 .4 1 24241 E-l C 
-6 . 4 1 CO 1 30E-1 C 
-6.4076498E-1C 
-6.4053331 E- 1C 
-6 .4C306 18E-1 C 
-6.4CC8347E-1C 
-6.39865C6E-1C 
-6 . 3965083E- 1 C 
-6.3944069E-1C 
-6 . 3923452E- 1 C 
-6.39C3223E-1C 
-6.3883372E-1C 
-6. 3863889 E- 1C 
-6. 3844765 E- 1C 
-6.3825991E-1C 
-6 . 38C7559E- 1 C 
-6.3789460E-1 C 
-6 . 3771 686E- 1 C 
-6.3754230E-1C 
-6.3737083E-1C 
-6.372C239E-1 C 
-6.37C3690E-1 C 
-6.3687429E-1C 
-6 .367 1449E-1 C 
-6.3655745E-1 C 
-6. 36403C9E-1 C 
-6. 3625136 E - 1C 
-6.361C219E-1C 
-6 . 3595553E- 1 C 
-6 . 3 58 1 1 32 E- 1 C 
-6.3566951 E-l C 
-6.3553004E-1C 
-6 . 3539285 E- 1 C 
-6.3525791 E- 1C 
-6. 3512516 E- 1C 
-6.3499455E-1C 
-6. 3486604 E-1C 
-6.3473958E-1C 
-6. 3461513 E- 1C 
-6 .3449264E-1 C 
-6.34 37208E- 1 C 
-6.3425340E- 1 C 
-6.34 13656E-1 C 

SOUND VELOCITY 
VERSUS 



1.0955184 M/S 
1.0948421 M/S 
1.0941818 M/S 
1.0935371 M/S 
1.0929075 M/S 
1.0922924 M./S 
1.0916914 M/S 
1.0911041 M/S 
1. 09C5300 M/S 
1.0899689 M/S 
1.0894202 M/S 
1.0888835 M/S 
1.0883586 M/S 
1.0878451 M/S 
1.0873427 M/S 
1. 0868510 M/S 
1.0863697 M/S 
1.0858985 M/S 
1. 0854372 M/S 
1.0849855 M/S 
1.0845430 M/S 
1.0841096 M/S 
1.0836850 M/S 
1.0832689 M/S 
1.0828612 M/S 
1.0824615 M/S 
1.082C697 M/S 
1.0816856 M/S 
1.0813C9C M/S 
1.0809396 M/S 
1.0805773 M/S 
1.0802219 M/S 
1.0798733 M/S 
1.0795311 M/S 
1.0791954 M/S 
1.0786659 M/S 
1.0785426 M/S 
1.0782250 M/S 
1.0779133 M/S 
1. 0776C73 M/S 
1 . 0773C67 M/S 
1.077C115 M/S 
1.0767215 M/S 
1.0764366 M/S 
1.0761567 M/S 
1.0758818 M/S 
1.0756115 M/S 
1.0753460 M/S 
1 . 075C849 M/S 
1.0748283 M/S 
1.0745761 M/S 
1.0743280 M/S 
1.074C842 M/S 
1.0738443 M/S 
1.0736085 M/S 
1.0733765 M/S 
1.0731403 M/S 
1.0729237 M/S 
1. 0727C29 M/S 
•1.0724856 M/S 
1.0722717 M/S 
1.0720613 M/S 
1.0718541 M/S 
1.0716502 M/S 
1.0714495 M/S 
1.0712519 M/S 



c'hJK WATER 
FREQ. 



TIME, 2 MINUTES AND 35 SECONDS 



WITH BUBBLES 
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APPENDIX VIII 



ENVELOPE DETECTOR AND ECHO COUNTER 

This system starts a Beckman 7360 one megacycle counter at the 
start of the main bang and stops it at the beginning of a preselected 
echo. Figure (IX-2) is a circuit diagram of the system. 

The system operation is as follows. The envelope detector 
passes the envelope of the echo train but does not pass the high 
frequency signal. The "and" circuit passes the echo train only 
between the main bang and the preselected echo for which we are 
measuring travel time. The amplifier and clipper convert the echo 
train into a rectangular pulse train with steep leading and trailing 
edges. The echo counter counts echos and sends a stop pulse to the 
one megacycle counter when the echo counter reaches the preselected 
count. At this time flip flop one is reset by this stop signal. 

The flip flop output turns off the clipper-amplifier. The next 
bang starts the one megacycle counter and the flip flo which in 
turn turns the amplifier-clipper back on. The process then 
repeats itself. 

The system as shown is plagued by two sources of error. First, 
as frequency increases the rise time of the detected wave pulse 
decreases thus affecting the counter reading. Second, as attenuation 
increases the rise time of the detected wave pulse increases 
affecting the counter reading. The combination of these effects 
caused velocity variation of five percent from 15 KC to 100 KC over 
two round trips in the laboratory "non-diBpersive" water. 
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By using a more elaborate clipping network with high sensitivity 
perhaps the above effects could be minimized. The wave shape of 
the detected pulse train limited the gain that could be used in the 
clipper as shown in figure (IX~1). This problem might be solved by 
using several clippers with different sensitivity. The most sensi- 
tive clipper would be turned on only when the last echo to be counted 
is received. This would enable later echos to be received. 

Receiving the 16th echo may decrease the variation of 0.6 percent 
and increasing the gain of this last echo may decrease the vari- 
ation further. Figure (IX-3) is a possible block diagram of the 
proposed improved system. Time did not permit construction of 
the improved system. 




FIGURE VIII-1 
DETECTOR OUTPUT 
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FIGURE ^7TTT -3 
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APPENDIX IX 



TIME MARKER 

This system presents a display on a dual trace oscilloscope 
similar to a radar "A" scope. A selected echo is superimposed on 
a step voltage. The step can be moved to any point on the echo. 

This step stops a one megacycle counter which was started by the 
main bang. 

The heart of the system is a DuMont 404 B pulse generator 
with delay. Delays from two microseconds to 12.5 milliseconds can 
be obtained in four ranges. With the dual trace oscilloscope 
synchronized to the start of the delayed pulse any echo up to the 
11th echo can be displayed in detail. The pulse width can be finely 
controlled with a four inch diameter knob providing 100 microseconds 
change for a 360 degree turn. The rise and fall times of the pulse 
are 14 nanoseconds maximum. Pulse widths from .05 to 105 micro= 
seconds in three ranges are available. 

The Beckman 7360 one megacycle counter is started by the gating 
voltage output of the tone burst generator. This gating voltage 
output has an amplitude of 24 volts, and width determined by the 
main bang width. The counter is stopped by the trailing edge of the 
pulse generator ouput. The counter can be started and stopped by 
a minimum plus or minus one volt pulse. 

Overall error due to instabilities in the counter, pulse 
generator, and oscilloscope are indicated in figure (iX-l). A 
mean deviation of about .011 percent is obtained from Figure 
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(IX-1) when using the seventh echo for measurements and .0094 
percent when using the eleventh echo. When taking actual measure- 
ments the error can be reduced by visual integration of the counter 
readout. 
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APPENDIX X 



Extracts from N.E.L. Tower log together with temperature and 
oxygen data. 

10 February 1965 

Tides: High 1746 at 3.2 ft,; Lows 1119 at 0.2 ft., 2202 at 2.4 ft 
Wave Height: 3g ft.; Water Clarity 15 ft., Slicks, None. 

Wind: Speed 0-6 kts.; Direction from 100° to 140°. 

Depth 10 ft.. Time 1330 - 1614. 

Temp.: 13.2 + 0.1°C; Oxygen Content 96 + 3$ saturation. 

Depth 10-45 ft.; Time 1835 - 1855. 

Temp. 13.25 + 0.05°C; Oxygen Content 101 + 3 % saturation 
Depth 45-10 ft.; Time 1918 - 1958 

Temp. 13.3 + 0.1°C; Oxygen Content 9 8+4# saturation 
Depth 10 ft.; Time 2141 - 2353. 

Temp. 13.4 + 0.1°C; Oxygen Content 100.5 + 1.5$ saturation 
12 February 1965 

Water Clarity, 20 ft.; Wind, Sp ee d o _3 knots; Direction from 095° 
Depth 10-45 ft.; Time 1156 - 1218 

Temp. 13.2 + 0.1°C; Oxygen Content 99.5 + 3.5$ saturation 
Depth 45-10 ft.; Time 1224 - 1247. 

Temp. 13-15 + 0.05°C.; Oxygen Content 95 + 2$ saturation 
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15 February 1965 

Tides j High 0805 at 7.2 ft.; Low 1512 at-1.9 ft.; High 2119 at 5.0 ft. 

Wave Height 4 ft.; Water Clarity 20 ft.; Slicks - None 

Wind; Speed 2 kts.; Direction from 020° 

Depth 10-45 ft.; Time 0925 - 0942 

Temp. 12.85 + 0.05°C.; Oxygen Content 98 + 2$ saturation 
Depth 45-10 ft.,; Time 0950 - 1008 

Temp. Isothermal at 12.9°C.; Oxygen Content 94.5 + 3.5$ saturation 
Depth 25 ft.; Time 1608 - 1639 

Temp. 1305 + 0.5°C.; Oxygen Content 94 + 1$ saturation 

Oxygen data is not available after this time since the oxygen probe 

lead split and water entered the system. 

16 February 1965 

Tides; Low 0252 at 0.8 ft.; High 0849 at 7.0 ft.; Low 1549 at-1.7 ft. 
Wave Height 2g ft.; Water Clarity 12 ft.; Slicks , Some. 

Wind; Speed 5 knots; Direction from 095° • 
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APPENDIX XI 



Effect of Scatterer Position on Signal 
Experimentation was carried out to determine the effect on 
backscattered sound from a 1.5 inch styrofoam sphere when it was 
placed in a number of positions within the region between the 
transducer and reflector. 

The sphere was suspended in the space at the positions shown 
in Figure XI-4 and the set of accompaning photos were taken to 
show the backscatter. A partial analysis follows i 
Pulse length = 75.0 cm = 2.46 ft = 0.5 msec 
Cycles per pulse (at 49.05 KC) = 24.5 
The region is broken into 7 approximately equal intervals as 
shown in Figure XI-4. The distance from the transducer to the 
reflector is 76.4 cm or ^ pulse length. 

.. I SPPPPE AT 

I Po$/7~/d)A/ /. 

I \ © 1 1 1 1 -» 

r t a 3 4 5 P 



Figure XI-1 

The first backscatter from sphere will appear at t = 4/7(25) = 14.3 
cycles having travelled to position 1 and back to the transducer 
but it will be blanked by the transmitted pulse for 3/7(25) = 10.7 
cycles. Although the reflected pulse exists for a full 25 cycles 
it will be seen only after the source pulse has passed by the 
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sphere, that is it will be seen for only t = 4 / 7 ( 25 ) = 14»3 
cycles. This is not completely seen in photo 1 because of the 
extra wide blanking pulse applied to the main bang interval. 

There is then a break of 3/7(25) = 10.7 cycles = 0.214 msec before 
the reflection from the relfector can be seenj this checks with 
photo 1. The reflector echo also lasts for 25 cycles. To be more 
exact, however, any echo will persist for a time equivalent to 
the length of the source pulse plus the length of or thickness of 
the insonified scattererj that is the sphere echo will persist for 
75 + 2.54(15) = 78.1 cm = 26 cycles. A breakdown of photo 1 is 
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Position 2: 



At position 3 (3/7 along the axis), the situation is shown 
in Figure XI -3. 



AT 

POS/7-/OA/ 2. 



I 1 1 0 1 1 » 1 

T / 2 3 4 S R 



Figure XI-3 
Sphere at Position 2 

The scattering from the sphere will begin at 6/7(25) = 

21.4 cycles = 0.429 msec. It will persist for 26 cycles as 
before and will be followed by a "no sound" interval of l/7(25) 

= 3 • 57 cycles = 0.0715 msec; this checks with photo 2. 

Position 3 • 

In this position (4/7 along the axis), the backscattered 
sound begins at 8/7(25) = 28.6 cycles = 0.561 msec. There is no 
blanking by the source pulse since it terminated 0.0714 msec 
earlier. The scattered sound continues until 15/7(0.5) = 1.07 
msec but will be overlapped from 15/7(0.5) = 1.07 msec or for the 
last 0.07 msec of the scattered pulse. 

The analysis can be continued in a similar manner for positions 
4 and 5. 



187 



Photos 6 through 8 show the backscatter for corresponding 
positions of the sphere. 

Taking the average amplitude of the scatter from position 
3 as a reference it was found that for the data taken the total 
deviation from this value throughout the region was approximately 

1056 , 
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